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ABSTRACT 

A systematic survey to establish the true incidence rate of the Blazhko modulation 
among short-period, fundamental-mode, Galactic field RR Lyrae stars has been ac- 
complished. The Konkoly Blazhko Survey (KBS) was initiated in 2004. Since then 
more than 750 nights of observation have been devoted to this project. A sample of 30 
RRab stars was extensively observed, and light-curve modulation was detected in 14 
cases. The 47% occurrence rate of the modulation is much larger than any previous 
estimate. The significant increase of the detected incidence rate is mostly due to the 
discovery of small-amplitude modulation. Half of the Blazhko variables in our sam- 
ple show modulation with so small amplitude that definitely have been missed in the 
previous surveys. We have found that the modulation can be very unstable in some 
cases, e.g. RY Com showed regular modulation only during one part of the observa- 
tions while during two seasons it had stable light curve with abrupt, small changes 
in the pulsation amplitude. This type of light-curve variability is also hard to detect 
in other Survey's data. The larger frequency of the light-curve modulation of RRab 
stars makes it even more important to find the still lacking explanation of the Blazhko 
phenomenon. 

The validity of the [Fe/H](P, Lp^i) relation using the mean light curves of Blazhko 
variables is checked in our sample. We have found that the formula gives accurate 
result for small-modulation-amplitude Blazhko stars, and this is also the case for large- 
modulation-amplitude stars if the light curve has complete phase coverage. However, 
if the data of large-modulation-amplitude Blazhko stars are not extended enough (e.g. 
less than 500 data points from less than 15 nights), the formula may give false result 
due to the distorted shape of the mean light curve used. 

Key words: stars: horizontal branch - stars: variables: other ~ stars: oscillations 
(including pulsations) - methods: data analysis - techniques: photometric 



1 INTRODUCTION 

The light-curve modulation of RR Lyrae stars, the so called 
Blazhko effect, is a hundred-year-old puzzle of stellar pul- 
sation. The Fourier spectra of the light curves of Blazhko 
variables are characterized by modulation frequency series 
appearing in the vicinity of the pulsation frequency (/o) and 

* Based on observations collected with the automatic 60 cm tele- 
scope of Konkoly Observatory, Budapest, Svabhegy 
f E-mail: jurcsik@konkoly.hu 



its harmonic frequencies (fc/o, fc > 1). Based on the location 
and number of the independent modulation frequency com- 
ponents some attempts have been made to introduce differ- 
ent subtypes of the mo d ulation (vi, V2, BLl , BL2, BL2x2, 
etc. lAlcock et al.l (120031 ): iMoskalik fc Porettil (|2003l ')'l. these 
classification schemes, however, do not indeed help to dis- 
close the unknown triggering mechanism of the modulation, 
which is most probably common in each Blazhko star. Both 
observational and theoretical efforts are still needed to un- 
derstand the modulation of RR Lyrae stars. Though the 
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incidence rate of the modulation was estimated to be about 
25-30% based on visual, photographic and photoelectric ob- 
servations of Galactic f ield and globular cluster variables 
ISmithl Il98ll . I1995I: [sieidl.,198a l. in a summary paper on 



RR Lyrae stars IPrestonI (196J) wrote: "There is no sharp 



distinction between singly and multiply periodic variables. 
When observed with sufficient precision or for long enough 
intervals of time, all RR Lyrae may be or may becom e mul- 
tiply periodic" - quoting iBalazs-Detre fc Detrd (|l962l l. 

More recently, the large C CD surveys e.g. MACHO 



dAlcock et a l.l 12003!') and OGLE ("Moskali k fc Porettil l2003l : 
[Mizerski ,2003(: ISoszvnski et all [2003: C oUinge et al.l I2OO6I I 
found a similar or even smaller fraction of variables ex- 
hibiting light-curve modulation than determined by photo- 
graphic and photoelectric observations. The incidence rates 
of the modulation in the Large Magellanic Cloud (LMC) 
and in the Galactic bulge were estimated to be 12-15% and 
23-27%, respectively. These statistics on Blazhko variables 
may, however, be biased due to the limits of the photometric 
accuracy and/or the deficiencies in data sampling. Observa- 
tion in longer wavelength band, where the amplitude of the 
modalation signal is small, may also reduce the detection 
rate of the modulation. The incidence rates determined so 
far are only l ower limits for the occ urrence of the modulation 
according to ICoIhnge et al] (|2006h . The metal deficiency of 
the variables was suggested to explain the sm all fraction of 
Blazhko stars in the LMC, but lSmolej (120051 ) could not de- 
tect any significant difference between the metallicities of 
RR Lyrae stars in the LMC and in t he Galactic bulge. 

Hippar cos JeSA 1997,:i NSVS (|Wozniak et al.ll2004 ). 
and ASAS (|Poimanskill2005l ) data bases also provide a large 
number of RR Lyrae observations but these data are not 
suitable to estimate the true percentage of variables showing 
light-curve modul ation. Based o n the NSVS an d ASAS data, 
IWils et al.1 (|2006l l and Szczvai cl fc Fabrvckvl 1200.7.) found 
modulations in an unrealistically small fraction, about 5% 
of RRab stars, RR Lyrae variables pulsating in the funda- 
mental mode, also designated as RRO stars. Many of the 
well-known Blazhko variables have not even been detected 
to show light-curve modulation in the Hipparcos, NSVS and 
ASAS data. 

No systematic survey that aimed to determine the fre- 
quency of light-curve modulation of RRab stars in the Galac- 
tic field has been performed previously. These variables are 
bright enough to obtain accurate light curves even with 
small telescopes. Earlier photometric programs concerning 
field RRab stars focused on either the light-curve variations 
of known Blazhko variables or th e multicolour light curves 



of large samples of variables (e.g. SturchI 19661: Fitch et al. 
19661: iBookmever et al.ll 19771 : lLublll977l : ISchmidt fc Reiswig 



19931 ) in order to determine their physical parameters (tem- 
perature, metallicity, reddening, etc.). Though these latter 
studies covered a large fraction of the brighter Galactic 
field variables, the time span, data number, phase cover- 
age, and photometric accuracy of the measurements were 
not suitable to detect long-period and/or small-amplitude 
light-curve variations. 

The knowledge of the true incidence rate of Blazhko 
stars has a great importance in finding the correct expla- 
nation of the phenomenon, as it shows how common the 
circumstances that favour the appearance of the modula- 
tion are. Exploiting the potential of our full access to an 



automated photometric telescope, we initiated a systematic 
study of RRab stars to get an estimate of the incidence 
rate of the modulation. Though our survey concerns a much 
more limited sample of variables than were measured in the 
MACHO and OGLE projects, the photometric accuracy and 
the denseness of our data sampling allow us to detect the 
occurrence of light-curve modulation more reliably than in 
those surveys. The main results concerning the incidence 
rate of the modulation of the Konkoly Blazhko Survey are 
summarized in the present paper. 

Results for variables with unstable light curves observed 
in the KBS that do not have data suitable to discuss in 
separate papers (RY Com, BD Her and FK Vul) are also 
published here. 

Comparing the photometric and spectroscopic metal- 
licities of Blazhko RRab stars of the KBS sample we also 
draw some conclusions on the correc tness of the photometric 
metallicity (jjurcsik fc Kovacd[l996l ) derived from the mean 
light curves of Blazhko variables. 



2 OBSERVATIONS AND TARGET SELECTION 

The observations of the KBS began in 2004 using an auto- 
mated 60 cm telescope (Budapest, Svabhegy) equipped with 
a Wright Instruments 750x1100 CCD camera (FoV 17'x24') 
and standard BVRcIc filters. Data were corrected for at- 
mospheric extinction and were transformed to the standard 
BVRcIc photometric system. We obtained accurate, ex- 
tended, multicolour light curves of fundamental-mode RR 
Lyrae stars with pulsation periods shorter than half a day. 
About 4000 hours of measurements were gathered for 30 
variables on ~ 750 nights during the past five years. Each 
variable was observed at least as long as it could be undoubt- 
edly decided whether the light curve was stable on about 
20-100 days time base or it was modulated. The typical ac- 
curacy of the observations was about 5-20 mmag, depend- 
ing mainly on the brightness of the object and on weather 
conditions. The apparent maximum V brightness of the ob- 
served variables was in the 10.5-13.5 mag range. The accu- 
racy, extension, and denseness of the data enabled us to 
detect light-curve modulations with small maximum bright- 
ness variations (some hundredths of a magnitude), i.e., with 
Fourier amplitudes of the modulation frequency components 
in the mmag regime. Data are utilized twofold; accurate mul- 
ticolour light curves of monoperiodic RRab stars and de- 
tailed analyses of individual Blazhko variables have already 
been and are going to be published. 

When selecting stars to observe, special care was taken 
to choose variables with good comparison stars (close to the 
variable in projected distance, brightness, and colour) in or- 
der to gain the possible most accurate photometry with our 
instrument. Magnitude differences of the variable relative to 
a single comparison star were determined and analyzed for 
most of the stars. The constancy of the comparison stars 
were verified by measuring several relatively bright stars in 
the fields. Details for the constancy of the comparison stars 
are given in the papers discussing the light c urve variabil- 
ity of the individual objects (see e.g.. Fig. 7 in I Jurcsik et al.l 
l2008d ). The ensemble mean brightnesses of five, three and 
two neighbouring stars were calculated, and used as compar- 
ison stars' magnitudes for UZ UMa, SS Cnc and DM Cyg, re- 
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spectively. W e applied the Image Subtraction Method (ISIS; 
lAlardll2003l ) for the photometry of CZ Lac in order to elim- 
inate the additional light of its close companion. 

We have observed only short-period RRab stars (P < 
0.5 days) to have a sample with somewhat homogeneous 
physical parameters, as most of the physical properties of 
RRab stars relate, to some degree, to the period of the pul- 
sation. The confinement of our survey to short-period vari- 
ables also made the observation effective, it assured to reach 
a good coverage of the entire possible sample within a few 
years of observations. No star with declination below 5° was 
selected to be observed in order to obtain extended data 
sets. 

The philosophy of the target selection was to choose 
variables with poor quality photometric observations avail- 
able and to clear up some questionable cases. We endeav- 
oured to obtain the possible most unaffected sample in 
order to get an unbiased estimate of the incidence rate 
of the modulation. None of the observed stars was defi- 
nitely known to show Blazhko modulation previously. In 
three cases (RR Gem, MW Lyr, and DM Cyg) light-curve 
modulation was announced based on photographic or vi- 
sual observations, but these results were in confiict with 
further observations and/or the reanalysis of the original 
data could not confirm the formerly detec ted modulation 
(|S6dor fc JurcsikI l2005l : ISodor et al.f l2007bl ) . Our observa- 
tions revealed that these stars do indeed show light-curve 
modulation, but with significantly different modulation pe- 
riod, amplitude, and type (amplitude or phase modulation) 
than that was found previo usly. RY Com was suspe cted to 
show Blazhko modulation (jjurcsik fc Kovacs|[l996h . how- 
ever, the scatter aroun d minimum brightness shown by the 
available observations (|joneslll966l : iBookmever et al.lll977l ) 
might have arisen from photometric inaccuracy. Four of the 
ta rgets (TZ Aur, SS Cnc, R R Gem, and FH Vul) were used 
bv lJurcsik fc Kovacsl (|l996l ) to calibrate the metallicity ver- 
sus light-curve parameters' relation assuming that their light 
curves were stable. Our accurate measurements show that, 
in fact, two of these stars, SS Cnc and RR Gem are Blazhko 
variables. 

The 30 RRab stars for which extended photometric data 
were obtained in the Konkoly Survey comprise about 60% 
and 80% of the possible targets of the Northern Hemisphere 
summer and winter seasons' RRab stars matching our selec- 
tion critetia. 



3 THE KONKOLY BLAZHKO SURVEY 

About half of the 30 observed RRab stars (14 variables, 47%) 
were found to exhibit light-curve modulation. This is a much 
larger percentage than was given by any previous estimate. 
None of the discovered Blazhko variables could have been 
undoubtedly identified to show light-curve modulation from 
the Hipparcos, NSVS, ASAS or any other previous photom- 
etry. As examples, our observations are compared to the 
published light curves for two RR Lyrae stars discovered to 
show Blazhko modulation in our survey in Figs. [1] and (2] 

AQ Lyrae, one of the stars measured in the sur- 
vey, was previously s upposed to have stable light curve. 
ICastellani et al. 1 (|1998D observed AQ Lyr on six consecutive 
nights in 1995 and found no light-curve modulation. The 
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Figure 1. Comparison of the Konkoly and previous photoelectric 
and CCD observations of AQ Lyr. The Konkoly observations ar e 
shown by dots, the CCD obs ervations of | Castellani et al.l lll99Sl) . 
and the photoelectric data of lSturchI l|l966l l are plotted by trian- 
gles and circles, respectively. Based on the previous observations 
no light-curve modulation of AQ Lyr could be detected. 



earlier photoelectric observations of SturchI (Il966l) matched 
the CCD hght curve observed by ICastellani et al.1 (|l998l ) 
within the limits of the photometric accuracy. Our extended 
data have reveled, however, that in fact AQ Lyr shows large- 
amplitude modulation. Although the variation in its max- 
imum V brightness is ~ 0.3 mag, its light curve does not 
change significantly on a few days time base due to the 65 
days length of its Blazhko period. In Fig.[T]our V light curve 
of AQ Lyr is compared to the previous observations. 

In Fig.[2]the Konkoly photometry is compared with the 
ASAS-3 and NSVS light curves of UZ Vir. The Konkoly light 
curve shows about 0.3 mag variation in maximum bright- 
ness, and corresponding to the Pbi = 68 d period of the 
modulation many triplet components can be identified with 
0.015 cd~^ (cycle per day) separation from the pulsation 
components in the prewhitened spectrum. Despite its strong 
maximum brightness variation, UZ Vir has not been found 
to be a Blazhko variable either from the ASAS or the NSVS 
data. These data show scattered light curves and very high 
noise level of the residual spectra. Although the /o — fm El 
modulation side-frequency appears in the residual spectrum 
of the ASAS data, this is not the largest amplitude frequency 
and there is no other signal in the residual spectrum with 
the same separation from any of the other pulsation compo- 
nents. 

Table [1] summarizes the observations and the results of 
the 30 RRab stars measured in the Konkoly Blazhko Survey. 
The identification of the comparison stars are given in the 



^ The modulation frequency fm = 1/Pbi- The primary modula- 
tion features in the Fourier spectra of Blazhko stars are triplets 
at kfo ~ fm; kfo; kfo + fm frequencies. 
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Figure 2. Comparison of the Konkoly and previous CCD light curves of UZ Vir. The residual spectra of the data prewhitcned for the 
pulsation light variation are shown in the right-hand panels. The largest amplitude signals in the residual spectrum of the Konkoly data 
are modulation components at 0.015 cd~^ separation from the pulsation frequency components and ±1 cd~^ alias frequencies of these 
modulation components. Although one of the modulation components (/o — /m) appears in the residual spectrum of the AS AS data, 
too, this is not the highest amplitude signal and no other significant peak with the same separation from the pulsation components can 
be identified. The scarce NSVS data show a scattered light curve. 



second column. The third, fourth, and fifth columns list the 
filters, the number of nights, and the time interval of the 
observations, respectively. Previous photometric light-curve 
information on the targets are given in the 'light-curve his- 
tory' co lumn. Pulsation perio ds are taken either from the 
GCVS (|Kholopov et al.lll985l ) or from our works (6th, 7th 
columns in Table [T|. The references given in the last column 



give information on the access to the Konkoly photometric 
datalfl 

The basic parameters of the modulation of the 14 
Blazhko stars discovered in the Konkoly Blazhko Survey 
are listed in Table [5] The modulation periods and ampli- 
tudes were determined from our observations. The modula- 

^ All the publications of the Konkoly Blazhko group and the 
Konkoly photometric data can also be downloaded from the web- 
site: [http : / / www, konkoly. hu/24/publications| . 
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Table 1. Summary of the Konkoly RR Lyrae Survey. Blazhko and monoperiodic variables are typeset in boldface and italics, respectively. 



Star 


Comparison Star 


Filters 


No._ 


b 

time" 


light-curve history'^ 


-Po [d] 


Ref._£ 

Po data 


AQ Lyr 


GSC2.3 N21I000423 


BV{I)c 


85 


460 


stable 


0.357134 


1 


1 


V759 Cyg 


GSC2.3 N2IA000839 


v{r)c 


74 


152 


n/a 


0.360014 


1 


1 


SS Cnc 


ensemble mean of 3 stars 


BV(RI)c 


35 


79 


stable 


0.367337 


2 


2 


EZ Cep 


GSC 4521-00784 


BV(I)c 


26 


202 


n/a 


0.3790035 


3 


3 


BK Cas 


CSC 4025-01395 


V(I)c 


11 


337 


n/a 


0.3902700 


3 


3 


BR Tau 


NA45000305 


BV{I)c 


115 


843 


n/a 


0.3905928 


1 


1 


TZ Aur 


BD +41 1609 


BV{RI)c 


13 


29 


stable 


0.3916746 


4 


5 


ET Per 


GSC 3671-01241 


BV(RI)c 


12 


23 


n/a 


0.3940135 


3 


3 


RR Gem 


BD -1-31 1547 


bv\ri)c 


63 


111 


stable /contradictory 


0.3972893 


6 


6 


MW Lyr 


GSC2.3 N0223233663 


Bv\l)c 


177 


361 


contradictory 


0.3976742 


7 


7 


V378 Per 


GSC2.3 NCGO000977 


V{I)c 


15 


96 


n/a 


0.3987208 


9 


9 


XY And 


GSC2.3 NBXO000560 


V(I)r 


64 


467 


n/a 


0.398725 


1 


1 


FH Vul 


GSC2.3 N2P8000417 


BViDc 


8 


108 


stable 


0.405413 


10 


10 


CN Lyr 


GSC2.3 N24S000237 


BV{I)c 


8 


59 


stable 


0.4113823 


4 


10 


DM Cyg 


GSC2.3 N0330220980; N03302207371 


BV(l\c 


81 


446 


contradictory 


0.419863 


11 


11 


BK And 


GSC2.3 N078000076 


BV(I)c 


20 


99 


n/a 


0.4216093 


9 


9 


CZ Lac 


Image Subtraction Method 


BV(RI)c 


116 


465 


n/a 


0.432174 


1 


1 


GI Gem 


GSC2.3 N8N9000652 


BV{I)c 


22 


92 


n/a 


0.4332664 


12 


12 


FK Vul 


GSC2.3 N2P8000417 


BV{I)c 


14 


83 


n/a 


0.4340529 


4 


13 


SW CVn 


BD -^37 2310 


BV{I)c 


10 


58 


n/a 


0.441671 


12 


12 


BH Aur 


GSC 02397-00378 


V(RI)c 


12 


19 


n/a 


0.4560898 


4 


5 


UU Boo 


GSC2.3 N6AZ000508 


BV{RI)c 


16 


396 


n/a 


0.4569339 


9 


9 


UZ Vir 


GSC2.3 N5IZ000213 


BV{I)c 


70 


474 


n/a 


0.4593925 


1 


1 


CG Peg 


GSC2.3 N2MC000574 


BV{I)c 


11 


95 


stable 


0.4671382 


4 


10 


UZ UMa 


ensemble mean of 5 stars 


V 


30 


115 


n/a 


0.4668413 


8 


8 


RY Com 


GSC2.3 N5CI000130 


BV{I)c 


98 


768 


n/a^ 


0.468951 


13 


13 


SU Leo 


GSC2.3 N6WV000233 


V{I)c 


12 


123 


n/a 


0.4722633 


12 


12 


BD Her 


GSC2.3 N2BS000412 


BV{I)c 


16 


69 


n/a 


0.4739064 


4 


13 


RZ Cam 


GSC2.3 N7T2000280 


BV{I)c 


17 


75 


n/a 


0.4804514 


12 


12 


TW Lyn 


GSC 02971-00853 


BV{RI)c 


17 


26 


n/a 


0.4818600 


4 


5 



Number of nights of the observations. 
Time interval of the observations, 
n/a: not enough photometric data. 

Ref erenc es: 1) forthcoming pap ers of the Konkoly RR Lyr ae g roup; 2)|jurcsilL.sLaLl ll2006bl'l: 3'l[S6 dor ct al | (l2007a^: 4'llKholo Dov ct al.l 
lll985^: 5:1 ljurcsik et all ll2006ari: 61 lJurcsik et all (l2005ah: 71 ljurcsik et al.l bOOSd) : 81 ISodor et'al... (.200^1: 91 ljurcsik et al. (JOOS^ 

1 01 iKun et al.l ||2008| 1 :~1 Jurcsik et al.l l|2009bl l: 121 ljurcsik et al.l l|2008d l: 131 this paper. 

^ ljurcsik &: Kovacj lll99(: 1 noted possible Blazhko modulation. The previous observations l|joneslll96ilBookmever e al.lll977ll show, 
however, large scatter at minimum light that may bo attributed to photometric inaccuracy. Sec also Fig. El 



tion amplitudes listed in Table [2] give different measures of 
the strength of the modulation. The full amplitude of the 
maximum-light variation {A{V)max) and the Fourier am- 
plitude of the largest modulation peak in V band {A{V)) 
are hsted. To be comparable with the OGLE-II results, the 
last column gives the Ic amplitude of the highest modula- 
tion peak in the vicinity of the main pulsation frequency 
{A{Ic))- This frequency is not always identical with the 
largest-amplitude modulation-frequency component, as our 
sample shows that the highest-amplitude modulation com- 
ponent appears in the vicinity of the first or second harmonic 
components of the pulsation (2/o, and 3/o) in about the half 
of the Blazhko variables. 

The observations of the 16 mo noperiodic RR a b stars 
of the s ample wer e pub lished by Jurcsik et al] (|2006al . 
l2008al lbh ; iKun et all (|2008[ l , and ISodor et al.l (|2007al l .These 
stars were observed on 8-26 nights, spanning 19-443 days 
time intervals. Within the limits of our observational accu- 



racy the light curves of these stars were stable. However, we 
cannot exclude the possibility that these stars may also show 
short- or long-period "micro-modulation" with 1-2 mmag 
amplitude of the modulation side frequencies, which is be- 
low our detecta bility limit. A ccording to the compilation of 
Blazhko stars bv lSmithl(il995l l about half of the short-period 
Blazhko stars have modulation periods longer than 50 days. 
Therefore, it cannot be excluded either that some of the 
stars claimed here to be monoperiodic may have long modu- 
lation periods, longer than 3-4 times the time interval of the 
observations. Moreover, there is a third possibility for the 
non-detection of the modulation of these stars, namely the 
te mporal behaviour of the modulation as shown for RR Gem 
bv lSodor etlj] (|2007bh and for RY Com in the next section 
of this paper. These stars can be regarded as monoperiodic 
only with these restrictions in mind. 

Detailed analyses of RR Gem, SS Cnc, and DM Cyg, 
three Blazhko variables showing small-amplitude light-curve 
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Table 2. Summary of the modulation properties of the 14 
Blazhko stars discovered in the Konkoly Blazhko Survey. 



Star Pbi [d] Abi [mag] _^ 

A(F)™, A{V) A{Ic) 



0.022 
0.009 
0.009 
0.007 
0.004 
0.056 
0.031 
0.006 
0.016 
0.04 
0.036 

0.005 
0.08 



°' j4(V)max: Full amplitude of the maximum light variation 
in V band; AiV): Fourier amplitude of the largest amplitude 
modulation frequency component in V band; A{Ic)'- Fourier 
amplitude of the largest amplitude modulation frequency com- 
ponent in the vicinity of /o in Iq band. 

complex multiperiodicity of the modulation is detected. 
'-■ modulation properties in 2007. 

m odulation that was not detected prev iously were published 
bv lJurcsik et al.1 (|2005al . [2006bl . bOOgbl ) . The doubly periodi c 
modulation of UZ UMa was shown by ISodor et all (|2006h . 
while the first really detailed analysis of the fight curve and 
colour behaviour of a large-modulatio n-amplitude Blazhko 
variab le, MW Lyrae, was presented by I Jurcsik et al.l (|2008d . 
l2009al ). The observations of AQ Lyr, V759 Cyg, BR Tau, 
XY And, CZ Lac and UZ Vir are extended enough to per- 
form detailed analyses of their Blazhko behaviour, which 
will be the subjects of further separate papers. Results for 
the remaining three modulated stars, RY Com, FK Vul, and 
BD Her are presented in this paper. 

Only the half of the new Blazhko variables (XY And, 
BD Her, CZ Lac, AQ Lyr, MW Lyr, UZ Vir, and FK Vul) 
have modulation amplitudes large enough to be detectable 
with less accurate and/or less extended observations. Nev- 
ertheless, even though the variation in maximum V bright- 
ness is larger than 0.15 mag and the largest-amplitude 
modulation-frequency component has an amplitude larger 
than 0.02 mag for these stars, they were not detected to show 
Blazhko modulation in the Hipparcos, ASAS and NSVS 
data. 

Based on the maximum timings collected in the GEOS 
data bas^ and our observations, we have found that there is 
no indication for any strong, irregular period change in any 
of the mono-periodic RRab stars in the sample. UU Boo, 
SW CVn and probably EZ Cep show steady period increase, 
while the period of FH Vul is probably decreasing. The peri- 
ods of the other 12 mono-periodic RRab stars are stable with 
very small if any changes detected. On the contrary, strong, 

^ http:/ /rr-lyr.ast.obs-mip.fr/dbrr/dbrr-V1.0J).php 



irregular changes are evident in some of the Blazhko stars 
(e.g, AQ Lyr, RR Gem, BD Her, V759 Cyg). We thus con- 
clude that the pulsation periods of Blazhko variables tend to 
be less stable than the periods of unmodulated RRab stars. 
A similar result was derived from t he study of t he period 
changes of RR Lyrae stars in M3 by ISzeidll (|l965h . 

3.1 The special case of RY Com 

We followed the light-curve variation of RY Com in three 
consecutive seasons, in 2007, 2008 and 2009. 

The Konkoly hght curve of RY Com between JD 
2 454 205 and 2 454 283 (2007 seaso r i) is compared with the 
available photoelectric (|joneslll966l : iBookmeyer et al.lll977l ) 
and CCD (NSVS, ASAS) observations in Fig.H The scatter 
of these latter measurements is much larger than the ampli- 
tude of the detected light-curve variation in the Konkoly 
observations. The residual spectra of the four data sets 
show different characteristics, after prewhitening the data 
for the mean pulsation hght curves. The highest peaks in the 
residual spectrum of the Konkoly data are modulation fre- 
quency components at /o±/m (/m ~ 0.03 cd~^) with small, 
^ 0.01 mag amplitudes, whereas the residual spectra of the 
ASAS, NSVS and the old photoelectric data show peaks 
higher than 0.04 mag at different frequencies. Peaks close to 
integer frequencies dominate the spectra of the old photo- 
electric and the NSVS data, indicating long-term trends in 
the observations, while the residual spectrum of the 4-year 
long ASAS data shows multiple peaks in the close vicin- 
ity of the pulsation frequency components. This is a typical 
feature if the pulsation period is changing during the time 
base of the observations. In the ASAS data this is indeed 
the case, the pulsation period of the variable was 0.00001 
d longer during the first 2 years of the ASAS observations 
than during the second part of the data. Any light-curve 
modulation of RY Com was definitely hidden in the noise 
and other biases in all the previous observations. 

Hoping to refine the modulation properties of RY Com, 
we continued its observation in 2008. However, in this sea- 
son we could not find any sign of a regular modulation, in- 
stead, the maximum brightness of the pulsation light curve 
dropped by about 0.08 mag in a very short period (about 
a week) at around JD 2 454 530. Though the scatter of the 
phased fight curves before and after JD 2 454 530 was larger 
than typical for our observations especially around minimum 
light, we could not find any definite modulation period in 
the data. The star behaved similarly in 2009, a sudden jump 
was detected again in maximum brightness in this season, 
too, at around JD 2 454 905, and no clear modulation was 
evident. The observed light-curve variations were consistent 
in the B, V, Ic bands in each season. In Fig.Uthe fight curve 
history of RY Com between 2007 and 2009 is documented. 

Fig. [5] shows the O — C variation of RY Com based on 
archive data collected in the GEOS data base and on our 
recent observations. The plots indicate that the pulsation 
period of RY Com varies irregularly, showing both abrupt 
and continuous changes. During the Konkoly observations a 
period increase of 0.00007 d occurred. Unfortunately, it can- 
not be resolved from our data whether the period change oc- 
curred abruptly at around the drop of maximum brightness 
(JD 2 454 530) or it was continuous. 

RY Com is the second star in our sample with tem- 
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Figure 3. Comparison of the Konkoly 2007 light curve and previous photoelectric and CCD observations of RY Com (right-hand 
panels). The residual spectra of the data prewhitened for the pulsation period in the — 5cd~^ frequency range and in the ±0.2cd^^ 
vicinity of the main pulsation frequency arc shown in the middle and right-hand panels, respectively. The residual spectra of the Konkoly 
observations show small-amplitude (< 0.01 mag) modulation-frequency components corresponding to a modulation period of 32 days. 
The residual spectra of the other data are dominated by noise and large-amplitude signals due to long-term trends in the data and to 
pulsation period changes of the star. 



porarily occurring modulation. For RR Gem it was found 
that during the 70 years of the observations the modula- 
tion properties of the star changed significantly and that the 
modul ation was not detectable in the 1970-80s (|S6dor et al.l 
l2007bl ). The photographic data of RR Gem showed that the 
change of the maximum brightness of the mean light curve 



was most probably connected to an abrupt period change. A 
0.00006 d period increase was followed by a 0.1 mag bright- 
ening of the Bpg mean maximum magnitude of RR Gem, 
while the 0.08 mag drop of the mean maximum V bright- 
ness of RY Com was accompanied by a 0.00007 d pulsation 
period increase in 2008. 
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Figure 4. The Konkoly V light curves of RY Com in the 2007, 2008 and 2009 observing seasons are shown. Crosses denote data belonging 
to the larger pulsation-amplitude phase of the 32 d modulation detected in the 2007 season and data of the first parts of the 2008 and 
2009 observations. Plus signs denote data when the maximum brightness of the star was fainter. The bottom panels show the comparison 
— check stars magnitude differences (GSC N5CI000130 — GSC3.2 N5CI000137). These data are also phased with the pulsation period 
of RY Com to demonstrate that the maximum-brightness variation of RY Com docs not arise from any defect of the comparison star's 
magnitudes. 
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Figure 6. Relative V magnitudes of FK Vul and BD Her versus Julian Date and phased with the pulsation and modulation periods are 
plotted in the top, middle, and bottom panels, respectively. In the middle panels the mean light curves, determined as appropriate order 
Fourier fits with the pulsation frequency to the data, are also shown. 



3.2 FK Vul and BD Her 

FK Vul and BD Her were observed in the course of the KBS 
in 2008. Both stars were observed on more than 10 nights. 
According to these observations large-amplitude modula- 
tions of the hght curves are evident (see Fig. |6}. The ob- 
servations indicate that the modulation periods of FK Vul 
and BD Her are about 56 and 22 days, respectively. As we 
do not plan to continue to observe these stars, and the ob- 



tained data are not extended enough to allow us to perform 
detailed analyses of these Blazhko variables, their observa- 
tions are presented here. 



The photometric observations and maximum timings of 
RY Com, FK Vul, and BD Her are available in the online 
version of the journal as Supplementary Material (see Ta- 
bles [S] and [4] as samples of the electronic data). 
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Table 3. Konkoly CCD observations of RY Com, FK Vul, and 
BD Her. Standard BVIc magnitude and B — V, V — Ic colour 
diflfcronccs arc given relative to the GSC2.3 stars N5CI000130 (13 
05 16.84 +23 16 54.0), N2PT000030 (20 52 40.50 +22 26 25.4), 
and N2BS000413 (18 50 32.84 +16 27 28.0) for RY Com, FK Vul, 
and BD Her, respectively. The full table is available as Supporting 
Information with the online version of this article. 
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Table 4. Maximum timings derived from the Konkoly CCD V 
observations for RY Com, FK Vul, and BD Her. The full table is 
published as Supporting Material. 
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4 COMPARISON WITH THE MACHO, 
OGLE-I, AND OGLE-II STATISTICS 

We compare our results w ith the MACHO LMC 
(|Alcock et all l2003l). OGLE-I jMoskahk fc Porettil l2003h . 
and OGLE-II (jCollinge et al.ll2006[ ) Galactic bulge statistics, 
as only these publications give enough details on the vari- 
ables. These surveys are based on the observations of 6158, 
150, and 1888 fundamental-mode RR Lyrae stars, among 
them 944, 51, and 544 have periods shorter than half a 
day. The detected incidence rates of the modulation of these 
short-period variables are 16%, 27%, and 30%, while 47% of 
the Konkoly sample exhibit the Blazhko effect. 

We do not make any distinction between the differ- 
ent types of modulation according to the location of the 
detected modulation frequency components {ui, U2, BLl, 
BL2, BL2x2, etc.), which were used in the other surveys. 
Our previous results showed that the vi, BLl variables, 
which are characterized by only one modulation frequency 
component close to the main pulsation frequency, exhibit, 
in fact, Blazhko type modulation (BL2, showing equidis- 
tant triplet frequencies) with high asymmetry in the a mpli- 
tudes of the side lobe frequencies (| Jurcsik et al.|[2005bl ). We 
have also shown that the side lobe frequencies of an am- 
plitude and phase modulated harmo nic signal can naturall y 
have highly asymmet ric amplitudes (jSzeidl fc Jurcsil3l2009^ . 
ICollinee et"alT(|2006l 'l obtained a significantly higher ratio of 
BL2 stars to BLl stars than what had been found in the 
previous surveys. This was attributed mainly to the robust- 
ness and sensitivity of the method they applied to detect 
additional symmetrical frequency components. Their result 
also indicates that the detection of equidistant frequency 
components on both sides of the pulsation frequency de- 
pends on the accuracy of the data and the efficiency of the 
analysis performed. It seems that characterizing the modu- 



MACHO 
LMC 



OGLE-I 
Galactic bulge 




.000 .025 .050 .075 .100 .125 .150 .175 

OGLE-II 
Galactic bulge 



.000 .025 ,050 .075 .100 .125 .150 .175 
Mcduiatlcn amplitude [mag] 



.000 .025 .050 .075 ,100 ,125 ,150 ,175 
Modulation ampiilude [mag] 



Figure 7. Modulation-amplitude distribution of the Blazhko 
variabl es with pulsation-p eriods short er than 0.5 days for the MA - 
CHO jAlcock et al.ll2003|). OG LE-I llMoskalik fc Porettil [20o3) . 
OGLE-lT JCoIlinge et al.lbooel ). and Konkoly surveys. There are 
148, 14, 166, and 14 short-period Blazhko stars detected in the 
four surveys. The modulation amplitude corresponds to the V 
Fourier amplitude of the largest-amplitude modulation-frequency 
component. The OGLE Iq amplitudes are multiplied by 1.58 in 
order to be comparable with the V band amplitudes. The MA- 
CHO and OGLE-I, II data show similar distributions, the modula- 
tion is the most frequent with 0.050-0.075 mag amplitudes accord- 
ing to these surveys. On the contrary, the more densely sampled 
and more accurate Konkoly data show that the modulation has 
small amplitude {A^^d < 0.025 mag) the most frequently. The 
Konkoly results warn that a large fraction of Blazhko staxs, those 
with small modulation amplitudes, has escaped detection previ- 
ously. 



lation with frequency doublets instead of triplets does not 
refiect real differences of the modulation. The 1/2, BL2x2 
variables show multiperiodic modulation, which is probably 
also a common feature of Blazhko stars. 

The distribution of the modulation amplitude of short- 
period (Ppuis < 0.5 days) Blazhko variables are shown for 
the MACHO, OGLE-I, OGLE-II, and Konkoly data in 
Fig. [7] The amplitude corresponds to the largest-amplitude 
modulation-frequency component in V band. In order to 
have homogeneous data, the Ic band amplitudes of the 
OGLE-I and OGLE-II data are transformed to V amplitudes 
accord ing to the formula: AiV) /A{Ic) = 1-58 (|jurcsik et al.l 
l2005d ). Here we stress again that amplitudes of modulation 
frequency components only in the vicinity of the main pul- 
sation frequency were published for the OGLE-II data, but 
the largest-amplitude modulation-frequency component oc- 
curs in the vicinity of the higher harmonics of the pulsation 
frequency in some Blazhko variables. 

While the MACHO and OGLE data show similar mod- 
ulation amplitude distributions with the most frequent am- 
plitude of the modulation being between 0.050-0.075 mag, 
small-amplitude modulations [AiV) < 0.025 mag) are the 
most frequent in the Konkoly sample. Blazhko variables with 
such small modulation amplitudes constitute only 1%, 7%, 
and 10% of the MACHO, OGLE-I, and OGLE-II Blazhko 
samples, respectively, but 50% of the modulated RRab stars 
discovered in the Konkoly Survey belong to this group. 



Blazhko modulation statistics of RRab stars 11 



The V and Ic Fourier amplitudes given in the last two 
columns in Table [2] can be directly compared with the mod- 
ulation amplitudes measured in the MACHO and OGLE-II 
surveys. In our sample of fourteen Blazhko stars, six vari- 
ables have A{V) ^ 0.015 mag and A{Ic) sj 0.009 mag mod- 
ulation amplitudes. Modulations with similarly small ampli- 
tudes were not detected in the MACHO and OGLE-I sur- 
veys at all, and were detected only in 2% of the short-period 
RRab stars of the OGLE-H data. 

If only Blazhko variables with at least one large- 
amplitude modulation-frequency component {A{V) > 
0.025 mag) were considered in the Konkoly data, the inci- 
dence rate of the modulation would be 23%, similarly to the 
MACHO and OGLE results. Consequently, the differences in 
the observed incidence rates of the modulation between our 
and the MACHO and OGLE statistics do not arise from the 
small number of objects and/or any bias in the target selec- 
tion of our survey but from the fact, that the latter surveys 
could not detect small-amplitude light-curve modulations. 

As our survey concerns only 30 objects, which is statis- 
tically a small number, it is important to verify the signifi- 
cance of the res ults obtained on this sample. In the OGLE- 
II data used bv lCollinge et all l|2006l 'l there were 544 RRab 
stars with Ppuis < 0.5 days and among them 30.5%, 166 vari- 
ables were detected to show light-curve modulation. Assum- 
ing that the true incidence rate of the modulation among 
RRab stars equals to this ratio, the likelihood of detecting 
14 out of 30 variables to show light-curve modulation would 
be only 4.6%. Therefore, the small sample size of the KBS 
cannot be accounted for the large percentage of Blazhko 
variables detected. 

We have also checked whether the selection procedure 
that ICoUinge et al.l (|2006l ) applied to detect Blazhko vari- 
ables could or could not identify small-amplitude modu- 
lations. The last column in Table [2] lists the Ic ampli- 
tudes of the highest modulation-peaks in the vicinity of 
the main pulsation frequency of the Blazhko variables ob- 
served in the Konkoly Survey. To classify a variable to BLl, 
this peak has to be as large as 2.4 times the noise level. 
For a BL2 variable two peaks with symmetrical frequency- 
separations are required, each with an amplitude larger than 

1. 1 times the noi s e leve l according to the selection applied 
bv lCollinge et al.l(|2006l ). The average noise level of the light 
curves of the short-period Blazhko variables in the OGLE- 
II sample is 16.7 mmag ± 8.3mmag. According to the data 
in Table [2] the A{Ic) amplitudes of the small-modulation- 
amplitude variables are about the half of the average noise 
level of the OGLE-II data. 

The probability of detecting RR Gem- type, small- 
amplitude modulation in the OGLE-II survey was also 
tested using different subsets of our observations. The same 
number of data points (250) as the OGLE-II light curves 
contain were randomly selected from the Ig band light 
curve of RR Gem. Following the method of ICoUinge et al.l 
l|2006l ) we searched for modulation frequency components 
in the vicinity of the main pulsation component in the 
test data sets. As RR Gem has a modulation period of 

7.2 days, modulation frequency components were searched 
in a wider frequency ran ge (/o ± 0.2 cd~^) than it was done 
bv lCollinge et al.l (|2006l ). RR Gem was detected only twice 
from 1000 random selecti ons to be a Blazhko variable ac- 
cording to the criteria of ICoUinge et al] (|2006l ). Therefore, 



we conclude that week modulations like RR Gem has were 
definitely missed in the OGLE-II survey. 

Another reason of the underestimation of the number of 
Blazhko vari ables in the OGLE-II data lies in the detection 
criteria that ICoUinge et al.l (|2006l ) applied. They searched 
only for modulations with Pmod > 10 days period. Taking 
into account that in our survey 2 of the 14 new Blazhko 
stars have modulation periods shorter than 10 days, this 
selection limit also leads to miss some fraction of the Blazhko 
variables. 

Concerning the MACHO data, the standard deviations 
of the residuals after the removal of the pulsation and mod- 
ulation components from the light curves are in the 0.03- 
0.20 mag range. This very high noise level makes the de- 
tection of small-amplitude modulations very unlikely. Using 
the light-curve solution of RR Gem and the timings of the 
MACHO observations we generated artificial time series to 
test the detection probability of RR Gem-like modulation 
in the MACHO data. To mimic the noise properties of the 
MACHO data, Gaussian noise was added with 0.03, 0.06, 
and 0.09 mag rms statistics. Fourier analysis of the test data 
(prewhitened for the pulsation) showed that the higher mod- 
ulation peak in the vicinity of the main pulsation frequency 
had an amplitude larger than 2a of the spectrum only if 
the smallest (0.03 mag) rms noise was added. If the rms was 
0.06 mag, there was none among 500 test time series that 
had a modulation peak reaching 2a amplitude. For compar- 
ison, only 5% of the MACHO light curves have rms residual 
smaller than 0.06 mag. 

Our survey is confined only to short-pulsation-period 
variables, which tends to bias the sample in favour of metal- 
rich disk RR Lyraes. Therefore our conclusion that the true 
incidence rate of the modulation is significantly larger than 
it was estimated in the previous surveys is valid only for 
this population of RRab stars. As the MACHO and OGLE 
surveys do not show, however, large differences in the occur- 
rence rates of the modulation at different pulsation periods, 
it is a sound guess that the situation is the same for longer 
period variables, too. To confirm this supposition accurate 
photometric data are needed for longer period RRab stars, 
as well. 



5 THE PHOTOMETRIC METALLICITY OF 
BLAZHKO VARIABLES 

The metallicity of RRab stars can be derived from the pe- 
riod of the pulsation and the tp^i epoch independent phase 
difference of the V light curve with similar accuracy to that 
determined from low di s persion spect roscopic observations 
(IJurcsik fc KovacE|[l996l :lK ovac 3 l2005l ). However, we do not 
know whether this relation also holds for the mean light 
curves of Blazhko variables. If the phase coverage of both the 
pulsation and the modulation cycles is not complete enough 
then the mean light curves of large-modulation-amplitude 
Blazhko variables can be seriously distorted, giving rise to 
spurious value of the Fourier parame t ers. F or instance, us- 
ing the ASAS data of SS For lKovacj (|2005l ) found that the 
photometric metallicity deviates significantly from the spec- 
trosc opic value, however, th e recent, extended photometric 
data (|Kolenberg et al]|2008l ) define a mean light curve with 
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Table 5. Comparison of spectroscopic and photometric [Fe/H] values of Blazhko stars. Variables with small and large modulation 
amplitudes are typeset in italics and boldface, respectively. 



star 


No. of V data 


[Fe/H]sp (weij 


jht)" 


[Fc/H],p " 


[Fc/H]phot 






Suntzeff et al. ( 1994 ~) 


Lavden (1994) 






33 Cnc 


1400 


-0.16 (4.0) 


0.13 (2.5) 


-0.03 


-0.12 


RR Gem 


3000 


-0.14 (4.5) 


-0.13 (2.0) 


-0.14 


-0.16 


DM Cyg 


3000 


-0.16 (1.0) 


0.07 (2.0) 


-0.01 


-0.01 


RY Com 


1000 


-0.64 (0.5) 


-1.38 (1.0) 


-1.13 


-1.16; -1.16; -1.08; -1.23; -l.ie'^ 


AQ Lyr 


1400 


-0.19 (2.0) 






-0.05 


XY And 


2600 




-0.68 (1.0) 




-0.68 


CZ Lac 


8000 


-0.13 (2.0) 






-0.26 


FK Vul'* 


500 




-0.71 (1.0) 




-0.37 


BD Her<* 


400 


-0.13 (2.0) 






-1.05 



" Spectroscopic data transforme d to the metallicity scale of the photometric [Fe/H] and their weights. 
Weighted means as specified in I Jurcsik &: Kovacsl (Il996f) . 

Derived from the mean light curves of the 2007 and 2008, 2009 data before and after JD 2454530 and JD 2454905, respectively. 
Scarced data. 



Fourier parameters well matching the [Fe/H] (P, 1/331) rela- 
tion. 

The extended data we obtained for Blazhko variables 
enable us to check the accuracy of the photometric metallic- 
ities derived from the mean light curves of these variables. 
Among the 14 new Blazhko variables nine have spectroscopic 
abundance determinations. Table [5] lists these variables and 
summarizes their spectroscopic and photometric metallici- 
ties. In this sample there are four and five stars showing 
small- and large-amplitude modulations, respectively. The 
photometric metallicities of the small-modulation-amplitude 
variables calculated from well defined mean fight curves are 
in very good agreement with the spectroscopic values. We 
thus conclude that the photometric metallicity of small- 
modulation- amplitude Blazhko variables is correct, and it 
has the same accuracy as for unmodulated RRab stars if 
the Ught variation is sampled properly. 

The situation for Blazhko variables showing large- 
amplitude modulations is, however, not so obvious. Unfor- 
tunately, the spectroscopic information on these stars is less 
certain than for the small-modulation-amplitude variables. 
For those variables that have good data coverage (AQ Lyr, 
XY And and CZ Lac) the spectroscopic and photometric 
metallicit ies agre e reaso nably well. A similar result was de- 
rived by ISmoled (|2005h too, using published photometric 
and spectroscopic data of five large-modulation-amplitude 
Blazhko stars. For the two large-modulation-amplitude vari- 
ables in our sample, which have only 400-600 measurements 
(FK Vul and BD Her), however, the discrepancy between 
the photometric and spectroscopic metallicities is large. The 
mean Ught curves used to calculate the [Fe/H] phot for FK Vul 
and BD Her are also drawn in the middle panels of Fig. [6] It 
can bee seen that the descending branch of the light curve 
of FK Vul is undersampled, which may result in improper 
shape of the mean light curve. Also, the mean light curve of 
BD Her, especially around minimum light, has an unusual 
shape most probably due to unequal data sampling. It is an 
open question, however, that if the light curves of FK Vul 
and BD Her were better sampled, then their photometric 
[Fe/H] would or would not match their spectroscopic val- 
ues better. The good agreement between the photometric 



and spectroscopic metallicities of those large-modulation- 
amplitude Blazhko variables which have good light-curve 
coverage indicates, however, that if the light variation is well 
sampled then the Fourier parameters of the mean light curve 
can be safely used to calculate the metallicity. Nevertheless, 
to make this statement better grounded, both extended pho- 
tometric data of other large-modulation-amplitude variables 
and spectroscopic metallicities of these variables (e.g. for 
MW Lyr) are needed. 



6 SUMMARY 

• We have shown that with the inclusion of small- 
amplitude modulations the incidence rate of the modula- 
tion increases to about 50% for short-period RRab stars. 
The Konkoly Blazhko Survey showed that, in fact, small- 
amplitude modulations are equally frequent as large ampli- 
tude ones. The distribution of the amplitudes of the modula- 
tion show an increase towards smaller amplitudes according 
to our data. It suggests that the light curves of those RRab 
stars that are now supposed to be stable may also be mod- 
ulated but with very small amplitude, below our present 
detection limit. Moreover, we may have missed detecting 
small-amplitude modulations with long modulation periods 
in our sample. Therefore, it is likely that the incidence rate 
of the modulation we detected is still a lower limit to the 
true fraction of the RRab stars showing the Blazhko effect. 
Modulations with amplitudes of the maximum brightness 
variation in the millimagnitude regime may still have es- 
caped detection. The observations of recent space missions 
(MOST, CoRoT, Kepler) will clear this question soon. 

• There is another embarrassing feature of the Blazhko 
modulation, namely its temporal occurrence as documented 
recentl y in the study of the long term behaviour of RR Gemi- 
norum (|S6dor et al.ll2007bl ') and in the case of RY Com in the 
present paper. The modulation of RR Lyrae i s also know n 
to cease for a short period in every fourth year (|Szeidllll988l ). 
Based on these results there is a good chance that the modu- 
lation is a common feature of the pulsation of fundamental- 
mode RR Lyrae stars, it may be even a universal property 
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of these variables. This possibility warns that understanding 
the Blazhko modulation is crucial in the knowledge of the 
pulsation of RR Lyrae stars. 

• Detecting small-amplitude light-curve modulation of 
RRa b stars i s also relevant from the theoretical point of 
view. iKovacd (|l995l ) used the lack of small-amplitude mod- 
ulations as an argument against the obli que magnetic rota - 
tor model of the Blazhko phenomenon ((Shibahashill200d ). 
He says "Assuming random distribution [of the aspect an- 
gle], we should see some small- and some large-amplitude 
modulations. This is apparently not the case." The discovery 
of small- amplitude modulation contradicts this argument. 
Whatever the cause of the Blazhko modulation, it manifests 
itself in the full possible range of amplitudes. At the same 
time, simultaneous multi-periodic modulation with similar 
modulation amplitudes were also detected in the KBS. The 
significant difference between the modulation periods ob- 
served in these stars is a strong argument against those 
explanations of the Blazhko phenomenon that connect the 
modulation period to the rotation of the star. 

• Based on our photometric data and spectroscopic in- 
formation on Blazhko stars we have found that the photo- 
metric metallicity determined from the mean light curves of 
Blazhko variables are correct if the light curve is properly 
sampled, i.e., the mean light curve is not distorted by uneven 
data sampling of the different Blazhko phases. 

In the Konkoly Blazhko Survey we obtained extended 
data on 30 RRab stars with fundamental periods shorter 
than 0.5 days, 14 of them showing the Blazhko effect. This is 
a significantly larger incidence rate of the modulation than 
found previously. We cannot exclude the possibility, how- 
ever, that the large percentage of Blazhko stars in the sample 
of the KBS is related to the short pulsation period of these 
stars. If this were the case, the incidence rate of the modu- 
lation would also depend on the metallicity of the stars as 
short period RRab stars are, in general, less metal poor than 
longer period variables. In order to decide whether there is a 
period and/or metallicity dependence of the occurrence rate 
of the modulation we initiated a similar survey focusing on 
RRab stars of longer periods in 2009. 
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ABSTRACT 

A systematic survey to establish the true incidence rate of the Blazhko modulation 
among short-period, fundamental-mode, Galactic field RR Lyrae stars has been ac- 
complished. The Konkoly Blazhko Survey (KBS) was initiated in 2004. Since then 
more than 750 nights of observation have been devoted to this project. A sample of 30 
RRab stars was extensively observed, and light-curve modulation was detected in 14 
cases. The 47% occurrence rate of the modulation is much larger than any previous 
estimate. The significant increase of the detected incidence rate is mostly due to the 
discovery of small-amplitude modulation. Half of the Blazhko variables in our sam- 
ple show modulation with so small amplitude that definitely have been missed in the 
previous surveys. We have found that the modulation can be very unstable in some 
cases, e.g. RY Com showed regular modulation only during one part of the observa- 
tions while during two seasons it had stable light curve with abrupt, small changes 
in the pulsation amplitude. This type of light-curve variability is also hard to detect 
in other Survey's data. The larger frequency of the light-curve modulation of RRab 
stars makes it even more important to find the still lacking explanation of the Blazhko 
phenomenon. 

The validity of the [Fe/H](P, Lp^i) relation using the mean light curves of Blazhko 
variables is checked in our sample. We have found that the formula gives accurate 
result for small-modulation-amplitude Blazhko stars, and this is also the case for large- 
modulation-amplitude stars if the light curve has complete phase coverage. However, 
if the data of large-modulation-amplitude Blazhko stars are not extended enough (e.g. 
less than 500 data points from less than 15 nights), the formula may give false result 
due to the distorted shape of the mean light curve used. 

Key words: stars: horizontal branch - stars: variables: other ~ stars: oscillations 
(including pulsations) - methods: data analysis - techniques: photometric 



1 INTRODUCTION 

The light-curve modulation of RR Lyrae stars, the so called 
Blazhko effect, is a hundred-year-old puzzle of stellar pul- 
sation. The Fourier spectra of the light curves of Blazhko 
variables are characterized by modulation frequency series 
appearing in the vicinity of the pulsation frequency (/o) and 

* Based on observations collected with the automatic 60 cm tele- 
scope of Konkoly Observatory, Budapest, Svabhegy 
f E-mail: jurcsik@konkoly.hu 



its harmonic frequencies (fc/o, fc > 1). Based on the location 
and number of the independent modulation frequency com- 
ponents some attempts have been made to introduce differ- 
ent subtypes of the mo d ulation (vi, V2, BLl , BL2, BL2x2, 
etc. lAlcock et al.l (120031 ): iMoskalik fc Porettil (|2003l ')'l. these 
classification schemes, however, do not indeed help to dis- 
close the unknown triggering mechanism of the modulation, 
which is most probably common in each Blazhko star. Both 
observational and theoretical efforts are still needed to un- 
derstand the modulation of RR Lyrae stars. Though the 
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incidence rate of the modulation was estimated to be about 
25-30% based on visual, photographic and photoelectric ob- 
servations of Galactic f ield and globular cluster variables 
ISmithl Il98ll . I1995I: [sieidl.,198a l. in a summary paper on 



RR Lyrae stars IPrestonI (196J) wrote: "There is no sharp 



distinction between singly and multiply periodic variables. 
When observed with sufficient precision or for long enough 
intervals of time, all RR Lyrae may be or may becom e mul- 
tiply periodic" - quoting iBalazs-Detre fc Detrd (|l962l l. 

More recently, the large C CD surveys e.g. MACHO 



dAlcock et a l.l 12003!') and OGLE ("Moskali k fc Porettil l2003l : 
[Mizerski ,2003(: ISoszvnski et all [2003: C oUinge et al.l I2OO6I I 
found a similar or even smaller fraction of variables ex- 
hibiting light-curve modulation than determined by photo- 
graphic and photoelectric observations. The incidence rates 
of the modulation in the Large Magellanic Cloud (LMC) 
and in the Galactic bulge were estimated to be 12-15% and 
23-27%, respectively. These statistics on Blazhko variables 
may, however, be biased due to the limits of the photometric 
accuracy and/or the deficiencies in data sampling. Observa- 
tion in longer wavelength band, where the amplitude of the 
modalation signal is small, may also reduce the detection 
rate of the modulation. The incidence rates determined so 
far are only l ower limits for the occ urrence of the modulation 
according to ICoIhnge et al] (|2006h . The metal deficiency of 
the variables was suggested to explain the sm all fraction of 
Blazhko stars in the LMC, but lSmolej (120051 ) could not de- 
tect any significant difference between the metallicities of 
RR Lyrae stars in the LMC and in t he Galactic bulge. 

Hippar cos JeSA 1997,:i NSVS (|Wozniak et al.ll2004 ). 
and ASAS (|Poimanskill2005l ) data bases also provide a large 
number of RR Lyrae observations but these data are not 
suitable to estimate the true percentage of variables showing 
light-curve modul ation. Based o n the NSVS an d ASAS data, 
IWils et al.1 (|2006l l and Szczvai cl fc Fabrvckvl 1200.7.) found 
modulations in an unrealistically small fraction, about 5% 
of RRab stars, RR Lyrae variables pulsating in the funda- 
mental mode, also designated as RRO stars. Many of the 
well-known Blazhko variables have not even been detected 
to show light-curve modulation in the Hipparcos, NSVS and 
ASAS data. 

No systematic survey that aimed to determine the fre- 
quency of light-curve modulation of RRab stars in the Galac- 
tic field has been performed previously. These variables are 
bright enough to obtain accurate light curves even with 
small telescopes. Earlier photometric programs concerning 
field RRab stars focused on either the light-curve variations 
of known Blazhko variables or th e multicolour light curves 



of large samples of variables (e.g. SturchI 19661: Fitch et al. 
19661: iBookmever et al.ll 19771 : lLublll977l : ISchmidt fc Reiswig 



19931 ) in order to determine their physical parameters (tem- 
perature, metallicity, reddening, etc.). Though these latter 
studies covered a large fraction of the brighter Galactic 
field variables, the time span, data number, phase cover- 
age, and photometric accuracy of the measurements were 
not suitable to detect long-period and/or small-amplitude 
light-curve variations. 

The knowledge of the true incidence rate of Blazhko 
stars has a great importance in finding the correct expla- 
nation of the phenomenon, as it shows how common the 
circumstances that favour the appearance of the modula- 
tion are. Exploiting the potential of our full access to an 



automated photometric telescope, we initiated a systematic 
study of RRab stars to get an estimate of the incidence 
rate of the modulation. Though our survey concerns a much 
more limited sample of variables than were measured in the 
MACHO and OGLE projects, the photometric accuracy and 
the denseness of our data sampling allow us to detect the 
occurrence of light-curve modulation more reliably than in 
those surveys. The main results concerning the incidence 
rate of the modulation of the Konkoly Blazhko Survey are 
summarized in the present paper. 

Results for variables with unstable light curves observed 
in the KBS that do not have data suitable to discuss in 
separate papers (RY Com, BD Her and FK Vul) are also 
published here. 

Comparing the photometric and spectroscopic metal- 
licities of Blazhko RRab stars of the KBS sample we also 
draw some conclusions on the correc tness of the photometric 
metallicity (jjurcsik fc Kovacd[l996l ) derived from the mean 
light curves of Blazhko variables. 



2 OBSERVATIONS AND TARGET SELECTION 

The observations of the KBS began in 2004 using an auto- 
mated 60 cm telescope (Budapest, Svabhegy) equipped with 
a Wright Instruments 750x1100 CCD camera (FoV 17'x24') 
and standard BVRcIc filters. Data were corrected for at- 
mospheric extinction and were transformed to the standard 
BVRcIc photometric system. We obtained accurate, ex- 
tended, multicolour light curves of fundamental-mode RR 
Lyrae stars with pulsation periods shorter than half a day. 
About 4000 hours of measurements were gathered for 30 
variables on ~ 750 nights during the past five years. Each 
variable was observed at least as long as it could be undoubt- 
edly decided whether the light curve was stable on about 
20-100 days time base or it was modulated. The typical ac- 
curacy of the observations was about 5-20 mmag, depend- 
ing mainly on the brightness of the object and on weather 
conditions. The apparent maximum V brightness of the ob- 
served variables was in the 10.5-13.5 mag range. The accu- 
racy, extension, and denseness of the data enabled us to 
detect light-curve modulations with small maximum bright- 
ness variations (some hundredths of a magnitude), i.e., with 
Fourier amplitudes of the modulation frequency components 
in the mmag regime. Data are utilized twofold; accurate mul- 
ticolour light curves of monoperiodic RRab stars and de- 
tailed analyses of individual Blazhko variables have already 
been and are going to be published. 

When selecting stars to observe, special care was taken 
to choose variables with good comparison stars (close to the 
variable in projected distance, brightness, and colour) in or- 
der to gain the possible most accurate photometry with our 
instrument. Magnitude differences of the variable relative to 
a single comparison star were determined and analyzed for 
most of the stars. The constancy of the comparison stars 
were verified by measuring several relatively bright stars in 
the fields. Details for the constancy of the comparison stars 
are given in the papers discussing the light c urve variabil- 
ity of the individual objects (see e.g.. Fig. 7 in I Jurcsik et al.l 
l2008d ). The ensemble mean brightnesses of five, three and 
two neighbouring stars were calculated, and used as compar- 
ison stars' magnitudes for UZ UMa, SS Cnc and DM Cyg, re- 
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spectively. W e applied the Image Subtraction Method (ISIS; 
lAlardll2003l ) for the photometry of CZ Lac in order to elim- 
inate the additional light of its close companion. 

We have observed only short-period RRab stars (P < 
0.5 days) to have a sample with somewhat homogeneous 
physical parameters, as most of the physical properties of 
RRab stars relate, to some degree, to the period of the pul- 
sation. The confinement of our survey to short-period vari- 
ables also made the observation effective, it assured to reach 
a good coverage of the entire possible sample within a few 
years of observations. No star with declination below 5° was 
selected to be observed in order to obtain extended data 
sets. 

The philosophy of the target selection was to choose 
variables with poor quality photometric observations avail- 
able and to clear up some questionable cases. We endeav- 
oured to obtain the possible most unaffected sample in 
order to get an unbiased estimate of the incidence rate 
of the modulation. None of the observed stars was defi- 
nitely known to show Blazhko modulation previously. In 
three cases (RR Gem, MW Lyr, and DM Cyg) light-curve 
modulation was announced based on photographic or vi- 
sual observations, but these results were in confiict with 
further observations and/or the reanalysis of the original 
data could not confirm the formerly detec ted modulation 
(|S6dor fc JurcsikI l2005l : ISodor et al.f l2007bl ) . Our observa- 
tions revealed that these stars do indeed show light-curve 
modulation, but with significantly different modulation pe- 
riod, amplitude, and type (amplitude or phase modulation) 
than that was found previo usly. RY Com was suspe cted to 
show Blazhko modulation (jjurcsik fc Kovacs|[l996h . how- 
ever, the scatter aroun d minimum brightness shown by the 
available observations (|joneslll966l : iBookmever et al.lll977l ) 
might have arisen from photometric inaccuracy. Four of the 
ta rgets (TZ Aur, SS Cnc, R R Gem, and FH Vul) were used 
bv lJurcsik fc Kovacsl (|l996l ) to calibrate the metallicity ver- 
sus light-curve parameters' relation assuming that their light 
curves were stable. Our accurate measurements show that, 
in fact, two of these stars, SS Cnc and RR Gem are Blazhko 
variables. 

The 30 RRab stars for which extended photometric data 
were obtained in the Konkoly Survey comprise about 60% 
and 80% of the possible targets of the Northern Hemisphere 
summer and winter seasons' RRab stars matching our selec- 
tion critetia. 



3 THE KONKOLY BLAZHKO SURVEY 

About half of the 30 observed RRab stars (14 variables, 47%) 
were found to exhibit light-curve modulation. This is a much 
larger percentage than was given by any previous estimate. 
None of the discovered Blazhko variables could have been 
undoubtedly identified to show light-curve modulation from 
the Hipparcos, NSVS, ASAS or any other previous photom- 
etry. As examples, our observations are compared to the 
published light curves for two RR Lyrae stars discovered to 
show Blazhko modulation in our survey in Figs. [1] and (2] 

AQ Lyrae, one of the stars measured in the sur- 
vey, was previously s upposed to have stable light curve. 
ICastellani et al. 1 (|1998D observed AQ Lyr on six consecutive 
nights in 1995 and found no light-curve modulation. The 
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Figure 1. Comparison of the Konkoly and previous photoelectric 
and CCD observations of AQ Lyr. The Konkoly observations ar e 
shown by dots, the CCD obs ervations of | Castellani et al.l lll99Sl) . 
and the photoelectric data of lSturchI l|l966l l are plotted by trian- 
gles and circles, respectively. Based on the previous observations 
no light-curve modulation of AQ Lyr could be detected. 



earlier photoelectric observations of SturchI (Il966l) matched 
the CCD hght curve observed by ICastellani et al.1 (|l998l ) 
within the limits of the photometric accuracy. Our extended 
data have reveled, however, that in fact AQ Lyr shows large- 
amplitude modulation. Although the variation in its max- 
imum V brightness is ~ 0.3 mag, its light curve does not 
change significantly on a few days time base due to the 65 
days length of its Blazhko period. In Fig.[T]our V light curve 
of AQ Lyr is compared to the previous observations. 

In Fig.[2]the Konkoly photometry is compared with the 
ASAS-3 and NSVS light curves of UZ Vir. The Konkoly light 
curve shows about 0.3 mag variation in maximum bright- 
ness, and corresponding to the Pbi = 68 d period of the 
modulation many triplet components can be identified with 
0.015 cd~^ (cycle per day) separation from the pulsation 
components in the prewhitened spectrum. Despite its strong 
maximum brightness variation, UZ Vir has not been found 
to be a Blazhko variable either from the ASAS or the NSVS 
data. These data show scattered light curves and very high 
noise level of the residual spectra. Although the /o — fm El 
modulation side-frequency appears in the residual spectrum 
of the ASAS data, this is not the largest amplitude frequency 
and there is no other signal in the residual spectrum with 
the same separation from any of the other pulsation compo- 
nents. 

Table [1] summarizes the observations and the results of 
the 30 RRab stars measured in the Konkoly Blazhko Survey. 
The identification of the comparison stars are given in the 



^ The modulation frequency fm = 1/Pbi- The primary modula- 
tion features in the Fourier spectra of Blazhko stars are triplets 
at kfo ~ fm; kfo; kfo + fm frequencies. 
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Figure 2. Comparison of the Konkoly and previous CCD light curves of UZ Vir. The residual spectra of the data prewhitcned for the 
pulsation light variation are shown in the right-hand panels. The largest amplitude signals in the residual spectrum of the Konkoly data 
are modulation components at 0.015 cd~^ separation from the pulsation frequency components and ±1 cd~^ alias frequencies of these 
modulation components. Although one of the modulation components (/o — /m) appears in the residual spectrum of the AS AS data, 
too, this is not the highest amplitude signal and no other significant peak with the same separation from the pulsation components can 
be identified. The scarce NSVS data show a scattered light curve. 



second column. The third, fourth, and fifth columns list the 
filters, the number of nights, and the time interval of the 
observations, respectively. Previous photometric light-curve 
information on the targets are given in the 'light-curve his- 
tory' co lumn. Pulsation perio ds are taken either from the 
GCVS (|Kholopov et al.lll985l ) or from our works (6th, 7th 
columns in Table [T|. The references given in the last column 



give information on the access to the Konkoly photometric 
datalfl 

The basic parameters of the modulation of the 14 
Blazhko stars discovered in the Konkoly Blazhko Survey 
are listed in Table [5] The modulation periods and ampli- 
tudes were determined from our observations. The modula- 

^ All the publications of the Konkoly Blazhko group and the 
Konkoly photometric data can also be downloaded from the web- 
site: [http : / / www, konkoly. hu/24/publications| . 
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Table 1. Summary of the Konkoly RR Lyrae Survey. Blazhko and monoperiodic variables are typeset in boldface and italics, respectively. 



Star 


Comparison Star 


Filters 


No._ 


b 

time" 


light-curve history'^ 


-Po [d] 


Ref._£ 

Po data 


AQ Lyr 


GSC2.3 N21I000423 


BV{I)c 


85 


460 


stable 


0.357134 


1 


1 


V759 Cyg 


GSC2.3 N2IA000839 


v{r)c 


74 


152 


n/a 


0.360014 


1 


1 


SS Cnc 


ensemble mean of 3 stars 


BV(RI)c 


35 


79 


stable 


0.367337 


2 


2 


EZ Cep 


GSC 4521-00784 


BV(I)c 


26 


202 


n/a 


0.3790035 


3 


3 


BK Cas 


CSC 4025-01395 


V(I)c 


11 


337 


n/a 


0.3902700 


3 


3 


BR Tau 


NA45000305 


BV{I)c 


115 


843 


n/a 


0.3905928 


1 


1 


TZ Aur 


BD +41 1609 


BV{RI)c 


13 


29 


stable 


0.3916746 


4 


5 


ET Per 


GSC 3671-01241 


BV(RI)c 


12 


23 


n/a 


0.3940135 


3 


3 


RR Gem 


BD -1-31 1547 


bv\ri)c 


63 


111 


stable /contradictory 


0.3972893 


6 


6 


MW Lyr 


GSC2.3 N0223233663 


Bv\l)c 


177 


361 


contradictory 


0.3976742 


7 


7 


V378 Per 


GSC2.3 NCGO000977 


V{I)c 


15 


96 


n/a 


0.3987208 


9 


9 


XY And 


GSC2.3 NBXO000560 


V(I)r 


64 


467 


n/a 


0.398725 


1 


1 


FH Vul 


GSC2.3 N2P8000417 


BViDc 


8 


108 


stable 


0.405413 


10 


10 


CN Lyr 


GSC2.3 N24S000237 


BV{I)c 


8 


59 


stable 


0.4113823 


4 


10 


DM Cyg 


GSC2.3 N0330220980; N03302207371 


BV(l\c 


81 


446 


contradictory 


0.419863 


11 


11 


BK And 


GSC2.3 N078000076 


BV(I)c 


20 


99 


n/a 


0.4216093 


9 


9 


CZ Lac 


Image Subtraction Method 


BV(RI)c 


116 


465 


n/a 


0.432174 


1 


1 


GI Gem 


GSC2.3 N8N9000652 


BV{I)c 


22 


92 


n/a 


0.4332664 


12 


12 


FK Vul 


GSC2.3 N2P8000417 


BV{I)c 


14 


83 


n/a 


0.4340529 


4 


13 


SW CVn 


BD -^37 2310 


BV{I)c 


10 


58 


n/a 


0.441671 


12 


12 


BH Aur 


GSC 02397-00378 


V(RI)c 


12 


19 


n/a 


0.4560898 


4 


5 


UU Boo 


GSC2.3 N6AZ000508 


BV{RI)c 


16 


396 


n/a 


0.4569339 


9 


9 


UZ Vir 


GSC2.3 N5IZ000213 


BV{I)c 


70 


474 


n/a 


0.4593925 


1 


1 


CG Peg 


GSC2.3 N2MC000574 


BV{I)c 


11 


95 


stable 


0.4671382 


4 


10 


UZ UMa 


ensemble mean of 5 stars 


V 


30 


115 


n/a 


0.4668413 


8 


8 


RY Com 


GSC2.3 N5CI000130 


BV{I)c 


98 


768 


n/a^ 


0.468951 


13 


13 


SU Leo 


GSC2.3 N6WV000233 


V{I)c 


12 


123 


n/a 


0.4722633 


12 


12 


BD Her 


GSC2.3 N2BS000412 


BV{I)c 


16 


69 


n/a 


0.4739064 


4 


13 


RZ Cam 


GSC2.3 N7T2000280 


BV{I)c 


17 


75 


n/a 


0.4804514 


12 


12 


TW Lyn 


GSC 02971-00853 


BV{RI)c 


17 


26 


n/a 


0.4818600 


4 


5 



Number of nights of the observations. 
Time interval of the observations, 
n/a: not enough photometric data. 

Ref erenc es: 1) forthcoming pap ers of the Konkoly RR Lyr ae g roup; 2)|jurcsilL.sLaLl ll2006bl'l: 3'l[S6 dor ct al | (l2007a^: 4'llKholo Dov ct al.l 
lll985^: 5:1 ljurcsik et all ll2006ari: 61 lJurcsik et all (l2005ah: 71 ljurcsik et al.l bOOSd) : 81 ISodor et'al... (.200^1: 91 ljurcsik et al. (JOOS^ 

1 01 iKun et al.l ||2008| 1 :~1 Jurcsik et al.l l|2009bl l: 121 ljurcsik et al.l l|2008d l: 131 this paper. 

^ ljurcsik &: Kovacj lll99(: 1 noted possible Blazhko modulation. The previous observations l|joneslll96ilBookmever e al.lll977ll show, 
however, large scatter at minimum light that may bo attributed to photometric inaccuracy. Sec also Fig. El 



tion amplitudes listed in Table [2] give different measures of 
the strength of the modulation. The full amplitude of the 
maximum-light variation {A{V)max) and the Fourier am- 
plitude of the largest modulation peak in V band {A{V)) 
are hsted. To be comparable with the OGLE-II results, the 
last column gives the Ic amplitude of the highest modula- 
tion peak in the vicinity of the main pulsation frequency 
{A{Ic))- This frequency is not always identical with the 
largest-amplitude modulation-frequency component, as our 
sample shows that the highest-amplitude modulation com- 
ponent appears in the vicinity of the first or second harmonic 
components of the pulsation (2/o, and 3/o) in about the half 
of the Blazhko variables. 

The observations of the 16 mo noperiodic RR a b stars 
of the s ample wer e pub lished by Jurcsik et al] (|2006al . 
l2008al lbh ; iKun et all (|2008[ l , and ISodor et al.l (|2007al l .These 
stars were observed on 8-26 nights, spanning 19-443 days 
time intervals. Within the limits of our observational accu- 



racy the light curves of these stars were stable. However, we 
cannot exclude the possibility that these stars may also show 
short- or long-period "micro-modulation" with 1-2 mmag 
amplitude of the modulation side frequencies, which is be- 
low our detecta bility limit. A ccording to the compilation of 
Blazhko stars bv lSmithl(il995l l about half of the short-period 
Blazhko stars have modulation periods longer than 50 days. 
Therefore, it cannot be excluded either that some of the 
stars claimed here to be monoperiodic may have long modu- 
lation periods, longer than 3-4 times the time interval of the 
observations. Moreover, there is a third possibility for the 
non-detection of the modulation of these stars, namely the 
te mporal behaviour of the modulation as shown for RR Gem 
bv lSodor etlj] (|2007bh and for RY Com in the next section 
of this paper. These stars can be regarded as monoperiodic 
only with these restrictions in mind. 

Detailed analyses of RR Gem, SS Cnc, and DM Cyg, 
three Blazhko variables showing small-amplitude light-curve 
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Table 2. Summary of the modulation properties of the 14 
Blazhko stars discovered in the Konkoly Blazhko Survey. 



Star Pbi [d] Abi [mag] _^ 

A(F)™, A{V) A{Ic) 



0.022 
0.009 
0.009 
0.007 
0.004 
0.056 
0.031 
0.006 
0.016 
0.04 
0.036 

0.005 
0.08 



°' j4(V)max: Full amplitude of the maximum light variation 
in V band; AiV): Fourier amplitude of the largest amplitude 
modulation frequency component in V band; A{Ic)'- Fourier 
amplitude of the largest amplitude modulation frequency com- 
ponent in the vicinity of /o in Iq band. 

complex multiperiodicity of the modulation is detected. 
'-■ modulation properties in 2007. 

m odulation that was not detected prev iously were published 
bv lJurcsik et al.1 (|2005al . [2006bl . bOOgbl ) . The doubly periodi c 
modulation of UZ UMa was shown by ISodor et all (|2006h . 
while the first really detailed analysis of the fight curve and 
colour behaviour of a large-modulatio n-amplitude Blazhko 
variab le, MW Lyrae, was presented by I Jurcsik et al.l (|2008d . 
l2009al ). The observations of AQ Lyr, V759 Cyg, BR Tau, 
XY And, CZ Lac and UZ Vir are extended enough to per- 
form detailed analyses of their Blazhko behaviour, which 
will be the subjects of further separate papers. Results for 
the remaining three modulated stars, RY Com, FK Vul, and 
BD Her are presented in this paper. 

Only the half of the new Blazhko variables (XY And, 
BD Her, CZ Lac, AQ Lyr, MW Lyr, UZ Vir, and FK Vul) 
have modulation amplitudes large enough to be detectable 
with less accurate and/or less extended observations. Nev- 
ertheless, even though the variation in maximum V bright- 
ness is larger than 0.15 mag and the largest-amplitude 
modulation-frequency component has an amplitude larger 
than 0.02 mag for these stars, they were not detected to show 
Blazhko modulation in the Hipparcos, ASAS and NSVS 
data. 

Based on the maximum timings collected in the GEOS 
data bas^ and our observations, we have found that there is 
no indication for any strong, irregular period change in any 
of the mono-periodic RRab stars in the sample. UU Boo, 
SW CVn and probably EZ Cep show steady period increase, 
while the period of FH Vul is probably decreasing. The peri- 
ods of the other 12 mono-periodic RRab stars are stable with 
very small if any changes detected. On the contrary, strong, 

^ http:/ /rr-lyr.ast.obs-mip.fr/dbrr/dbrr-V1.0J).php 



irregular changes are evident in some of the Blazhko stars 
(e.g, AQ Lyr, RR Gem, BD Her, V759 Cyg). We thus con- 
clude that the pulsation periods of Blazhko variables tend to 
be less stable than the periods of unmodulated RRab stars. 
A similar result was derived from t he study of t he period 
changes of RR Lyrae stars in M3 by ISzeidll (|l965h . 

3.1 The special case of RY Com 

We followed the light-curve variation of RY Com in three 
consecutive seasons, in 2007, 2008 and 2009. 

The Konkoly hght curve of RY Com between JD 
2 454 205 and 2 454 283 (2007 seaso r i) is compared with the 
available photoelectric (|joneslll966l : iBookmeyer et al.lll977l ) 
and CCD (NSVS, ASAS) observations in Fig.H The scatter 
of these latter measurements is much larger than the ampli- 
tude of the detected light-curve variation in the Konkoly 
observations. The residual spectra of the four data sets 
show different characteristics, after prewhitening the data 
for the mean pulsation hght curves. The highest peaks in the 
residual spectrum of the Konkoly data are modulation fre- 
quency components at /o±/m (/m ~ 0.03 cd~^) with small, 
^ 0.01 mag amplitudes, whereas the residual spectra of the 
ASAS, NSVS and the old photoelectric data show peaks 
higher than 0.04 mag at different frequencies. Peaks close to 
integer frequencies dominate the spectra of the old photo- 
electric and the NSVS data, indicating long-term trends in 
the observations, while the residual spectrum of the 4-year 
long ASAS data shows multiple peaks in the close vicin- 
ity of the pulsation frequency components. This is a typical 
feature if the pulsation period is changing during the time 
base of the observations. In the ASAS data this is indeed 
the case, the pulsation period of the variable was 0.00001 
d longer during the first 2 years of the ASAS observations 
than during the second part of the data. Any light-curve 
modulation of RY Com was definitely hidden in the noise 
and other biases in all the previous observations. 

Hoping to refine the modulation properties of RY Com, 
we continued its observation in 2008. However, in this sea- 
son we could not find any sign of a regular modulation, in- 
stead, the maximum brightness of the pulsation light curve 
dropped by about 0.08 mag in a very short period (about 
a week) at around JD 2 454 530. Though the scatter of the 
phased fight curves before and after JD 2 454 530 was larger 
than typical for our observations especially around minimum 
light, we could not find any definite modulation period in 
the data. The star behaved similarly in 2009, a sudden jump 
was detected again in maximum brightness in this season, 
too, at around JD 2 454 905, and no clear modulation was 
evident. The observed light-curve variations were consistent 
in the B, V, Ic bands in each season. In Fig.Uthe fight curve 
history of RY Com between 2007 and 2009 is documented. 

Fig. [5] shows the O — C variation of RY Com based on 
archive data collected in the GEOS data base and on our 
recent observations. The plots indicate that the pulsation 
period of RY Com varies irregularly, showing both abrupt 
and continuous changes. During the Konkoly observations a 
period increase of 0.00007 d occurred. Unfortunately, it can- 
not be resolved from our data whether the period change oc- 
curred abruptly at around the drop of maximum brightness 
(JD 2 454 530) or it was continuous. 

RY Com is the second star in our sample with tem- 



AQ Lyr 


64.9^ 
16.0^' 


0.30 


0.032 


V759 Cyg 


0.12 


0.014 


SS Cnc 


5.3 


0.09 


0.015 


BR Tau 


19.3 


0.13 


0.012 


RR Gem 


7.2 


0.09 


0.007 


MW Lyr 


16.5 


0.45 


0.090 


XY And 


41.4 


0.20 


0.050 


DM Cyg 


10.6 


0.07 


0.010 


CZ Lac 


14.6/18.6 


0.46 


0.027 


FK Vul 


~56 


0.3 


0.06 


UZ Vir 


68.2 


0.36 


0.053 


UZ UMa 


26.7/143: 


0.14 


0.020 


RY Com^ 


32: 


0.06 


0.008 


BD Her" 


~22 


0.5 


0.12 
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Figure 3. Comparison of the Konkoly 2007 light curve and previous photoelectric and CCD observations of RY Com (right-hand 
panels). The residual spectra of the data prewhitened for the pulsation period in the — 5cd~^ frequency range and in the ±0.2cd^^ 
vicinity of the main pulsation frequency arc shown in the middle and right-hand panels, respectively. The residual spectra of the Konkoly 
observations show small-amplitude (< 0.01 mag) modulation-frequency components corresponding to a modulation period of 32 days. 
The residual spectra of the other data are dominated by noise and large-amplitude signals due to long-term trends in the data and to 
pulsation period changes of the star. 



porarily occurring modulation. For RR Gem it was found 
that during the 70 years of the observations the modula- 
tion properties of the star changed significantly and that the 
modul ation was not detectable in the 1970-80s (|S6dor et al.l 
l2007bl ). The photographic data of RR Gem showed that the 
change of the maximum brightness of the mean light curve 



was most probably connected to an abrupt period change. A 
0.00006 d period increase was followed by a 0.1 mag bright- 
ening of the Bpg mean maximum magnitude of RR Gem, 
while the 0.08 mag drop of the mean maximum V bright- 
ness of RY Com was accompanied by a 0.00007 d pulsation 
period increase in 2008. 



8 J. Jurcsik et al. 



2007 2008 2009 





-1.4 


1 1 1 1 


i 1 i 1 1 1 


1 1 1 1 i 1 1 1 i 1 i 1 1 1 1 1 1 i 
" : s I" 

■ ■ ■: . t f , ti . . . ; t 


i 1 1 1 1 1 i 1 1 1 

■ ^ : =■ 


;om-C1) 


-1.2 
-1 


* J i ; 




■ il: I ■ 1 ;!!;!;: H ■ 




ta V (RY C 


-0.8 
-0.6 


* 4 - ^ : 
- 1-11 \ 




.'',14 ' * . * * • • 


" . "< 1 i*. . -j J * : " 

S ; 'I •• .-'tit ' 


(U 

■a 


-0.4 
-0.2 


111' 
1 . . . 


* * / + « s 
.• . , " 
|; . ■ i - 

.Mi i : ! ■ 

1 . . . . 


■ ! i| M i i i i 

. 1 .... 1 .... 1 .... 1 . 


■ 1 i% t |i« h ■ 
~. , 1 , , , . 1 , , ~ 



54200 54250 54500 54550 54600 54650 54900 54950 

JD-2400000 JD-2400000 JD-2400000 




0.5 1 0.5 1 0.5 1 



pulsation phase pulsation phase pulsation phase 



Figure 4. The Konkoly V light curves of RY Com in the 2007, 2008 and 2009 observing seasons are shown. Crosses denote data belonging 
to the larger pulsation-amplitude phase of the 32 d modulation detected in the 2007 season and data of the first parts of the 2008 and 
2009 observations. Plus signs denote data when the maximum brightness of the star was fainter. The bottom panels show the comparison 
— check stars magnitude differences (GSC N5CI000130 — GSC3.2 N5CI000137). These data are also phased with the pulsation period 
of RY Com to demonstrate that the maximum-brightness variation of RY Com docs not arise from any defect of the comparison star's 
magnitudes. 
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recent Konkoly data. The right-hand panel is an enlarged part of the O — C diagram for the last three years using the Konkoly data. 
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Figure 6. Relative V magnitudes of FK Vul and BD Her versus Julian Date and phased with the pulsation and modulation periods are 
plotted in the top, middle, and bottom panels, respectively. In the middle panels the mean light curves, determined as appropriate order 
Fourier fits with the pulsation frequency to the data, are also shown. 



3.2 FK Vul and BD Her 

FK Vul and BD Her were observed in the course of the KBS 
in 2008. Both stars were observed on more than 10 nights. 
According to these observations large-amplitude modula- 
tions of the hght curves are evident (see Fig. |6}. The ob- 
servations indicate that the modulation periods of FK Vul 
and BD Her are about 56 and 22 days, respectively. As we 
do not plan to continue to observe these stars, and the ob- 



tained data are not extended enough to allow us to perform 
detailed analyses of these Blazhko variables, their observa- 
tions are presented here. 



The photometric observations and maximum timings of 
RY Com, FK Vul, and BD Her are available in the online 
version of the journal as Supplementary Material (see Ta- 
bles [S] and [4] as samples of the electronic data). 
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Table 3. Konkoly CCD observations of RY Com, FK Vul, and 
BD Her. Standard BVIc magnitude and B — V, V — Ic colour 
diflfcronccs arc given relative to the GSC2.3 stars N5CI000130 (13 
05 16.84 +23 16 54.0), N2PT000030 (20 52 40.50 +22 26 25.4), 
and N2BS000413 (18 50 32.84 +16 27 28.0) for RY Com, FK Vul, 
and BD Her, respectively. The full table is available as Supporting 
Information with the online version of this article. 



HJD 

2400000+ 


delta mag 


filter / colour 


star 


54205.28742 


-1.911 


B 


RY.Com 


54205.29319 


-1.862 


B 


RY.Com 


54205.29978 


-1.835 


B 


RY.Com 



Table 4. Maximum timings derived from the Konkoly CCD V 
observations for RY Com, FK Vul, and BD Her. The full table is 
published as Supporting Material. 



maximum time 


star 


54211.368 


RY.Com 


54218.404 


RY.Com 


54227.315 


RY.Com 



4 COMPARISON WITH THE MACHO, 
OGLE-I, AND OGLE-II STATISTICS 

We compare our results w ith the MACHO LMC 
(|Alcock et all l2003l). OGLE-I jMoskahk fc Porettil l2003h . 
and OGLE-II (jCollinge et al.ll2006[ ) Galactic bulge statistics, 
as only these publications give enough details on the vari- 
ables. These surveys are based on the observations of 6158, 
150, and 1888 fundamental-mode RR Lyrae stars, among 
them 944, 51, and 544 have periods shorter than half a 
day. The detected incidence rates of the modulation of these 
short-period variables are 16%, 27%, and 30%, while 47% of 
the Konkoly sample exhibit the Blazhko effect. 

We do not make any distinction between the differ- 
ent types of modulation according to the location of the 
detected modulation frequency components {ui, U2, BLl, 
BL2, BL2x2, etc.), which were used in the other surveys. 
Our previous results showed that the vi, BLl variables, 
which are characterized by only one modulation frequency 
component close to the main pulsation frequency, exhibit, 
in fact, Blazhko type modulation (BL2, showing equidis- 
tant triplet frequencies) with high asymmetry in the a mpli- 
tudes of the side lobe frequencies (| Jurcsik et al.|[2005bl ). We 
have also shown that the side lobe frequencies of an am- 
plitude and phase modulated harmo nic signal can naturall y 
have highly asymmet ric amplitudes (jSzeidl fc Jurcsil3l2009^ . 
ICollinee et"alT(|2006l 'l obtained a significantly higher ratio of 
BL2 stars to BLl stars than what had been found in the 
previous surveys. This was attributed mainly to the robust- 
ness and sensitivity of the method they applied to detect 
additional symmetrical frequency components. Their result 
also indicates that the detection of equidistant frequency 
components on both sides of the pulsation frequency de- 
pends on the accuracy of the data and the efficiency of the 
analysis performed. It seems that characterizing the modu- 
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Figure 7. Modulation-amplitude distribution of the Blazhko 
variabl es with pulsation-p eriods short er than 0.5 days for the MA - 
CHO jAlcock et al.ll2003|). OG LE-I llMoskalik fc Porettil [20o3) . 
OGLE-lT JCoIlinge et al.lbooel ). and Konkoly surveys. There are 
148, 14, 166, and 14 short-period Blazhko stars detected in the 
four surveys. The modulation amplitude corresponds to the V 
Fourier amplitude of the largest-amplitude modulation-frequency 
component. The OGLE Iq amplitudes are multiplied by 1.58 in 
order to be comparable with the V band amplitudes. The MA- 
CHO and OGLE-I, II data show similar distributions, the modula- 
tion is the most frequent with 0.050-0.075 mag amplitudes accord- 
ing to these surveys. On the contrary, the more densely sampled 
and more accurate Konkoly data show that the modulation has 
small amplitude {A^^d < 0.025 mag) the most frequently. The 
Konkoly results warn that a large fraction of Blazhko staxs, those 
with small modulation amplitudes, has escaped detection previ- 
ously. 



lation with frequency doublets instead of triplets does not 
refiect real differences of the modulation. The 1/2, BL2x2 
variables show multiperiodic modulation, which is probably 
also a common feature of Blazhko stars. 

The distribution of the modulation amplitude of short- 
period (Ppuis < 0.5 days) Blazhko variables are shown for 
the MACHO, OGLE-I, OGLE-II, and Konkoly data in 
Fig. [7] The amplitude corresponds to the largest-amplitude 
modulation-frequency component in V band. In order to 
have homogeneous data, the Ic band amplitudes of the 
OGLE-I and OGLE-II data are transformed to V amplitudes 
accord ing to the formula: AiV) /A{Ic) = 1-58 (|jurcsik et al.l 
l2005d ). Here we stress again that amplitudes of modulation 
frequency components only in the vicinity of the main pul- 
sation frequency were published for the OGLE-II data, but 
the largest-amplitude modulation-frequency component oc- 
curs in the vicinity of the higher harmonics of the pulsation 
frequency in some Blazhko variables. 

While the MACHO and OGLE data show similar mod- 
ulation amplitude distributions with the most frequent am- 
plitude of the modulation being between 0.050-0.075 mag, 
small-amplitude modulations [AiV) < 0.025 mag) are the 
most frequent in the Konkoly sample. Blazhko variables with 
such small modulation amplitudes constitute only 1%, 7%, 
and 10% of the MACHO, OGLE-I, and OGLE-II Blazhko 
samples, respectively, but 50% of the modulated RRab stars 
discovered in the Konkoly Survey belong to this group. 
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The V and Ic Fourier amplitudes given in the last two 
columns in Table [2] can be directly compared with the mod- 
ulation amplitudes measured in the MACHO and OGLE-II 
surveys. In our sample of fourteen Blazhko stars, six vari- 
ables have A{V) ^ 0.015 mag and A{Ic) sj 0.009 mag mod- 
ulation amplitudes. Modulations with similarly small ampli- 
tudes were not detected in the MACHO and OGLE-I sur- 
veys at all, and were detected only in 2% of the short-period 
RRab stars of the OGLE-H data. 

If only Blazhko variables with at least one large- 
amplitude modulation-frequency component {A{V) > 
0.025 mag) were considered in the Konkoly data, the inci- 
dence rate of the modulation would be 23%, similarly to the 
MACHO and OGLE results. Consequently, the differences in 
the observed incidence rates of the modulation between our 
and the MACHO and OGLE statistics do not arise from the 
small number of objects and/or any bias in the target selec- 
tion of our survey but from the fact, that the latter surveys 
could not detect small-amplitude light-curve modulations. 

As our survey concerns only 30 objects, which is statis- 
tically a small number, it is important to verify the signifi- 
cance of the res ults obtained on this sample. In the OGLE- 
II data used bv lCollinge et all l|2006l 'l there were 544 RRab 
stars with Ppuis < 0.5 days and among them 30.5%, 166 vari- 
ables were detected to show light-curve modulation. Assum- 
ing that the true incidence rate of the modulation among 
RRab stars equals to this ratio, the likelihood of detecting 
14 out of 30 variables to show light-curve modulation would 
be only 4.6%. Therefore, the small sample size of the KBS 
cannot be accounted for the large percentage of Blazhko 
variables detected. 

We have also checked whether the selection procedure 
that ICoUinge et al.l (|2006l ) applied to detect Blazhko vari- 
ables could or could not identify small-amplitude modu- 
lations. The last column in Table [2] lists the Ic ampli- 
tudes of the highest modulation-peaks in the vicinity of 
the main pulsation frequency of the Blazhko variables ob- 
served in the Konkoly Survey. To classify a variable to BLl, 
this peak has to be as large as 2.4 times the noise level. 
For a BL2 variable two peaks with symmetrical frequency- 
separations are required, each with an amplitude larger than 

1. 1 times the noi s e leve l according to the selection applied 
bv lCollinge et al.l(|2006l ). The average noise level of the light 
curves of the short-period Blazhko variables in the OGLE- 
II sample is 16.7 mmag ± 8.3mmag. According to the data 
in Table [2] the A{Ic) amplitudes of the small-modulation- 
amplitude variables are about the half of the average noise 
level of the OGLE-II data. 

The probability of detecting RR Gem- type, small- 
amplitude modulation in the OGLE-II survey was also 
tested using different subsets of our observations. The same 
number of data points (250) as the OGLE-II light curves 
contain were randomly selected from the Ig band light 
curve of RR Gem. Following the method of ICoUinge et al.l 
l|2006l ) we searched for modulation frequency components 
in the vicinity of the main pulsation component in the 
test data sets. As RR Gem has a modulation period of 

7.2 days, modulation frequency components were searched 
in a wider frequency ran ge (/o ± 0.2 cd~^) than it was done 
bv lCollinge et al.l (|2006l ). RR Gem was detected only twice 
from 1000 random selecti ons to be a Blazhko variable ac- 
cording to the criteria of ICoUinge et al] (|2006l ). Therefore, 



we conclude that week modulations like RR Gem has were 
definitely missed in the OGLE-II survey. 

Another reason of the underestimation of the number of 
Blazhko vari ables in the OGLE-II data lies in the detection 
criteria that ICoUinge et al.l (|2006l ) applied. They searched 
only for modulations with Pmod > 10 days period. Taking 
into account that in our survey 2 of the 14 new Blazhko 
stars have modulation periods shorter than 10 days, this 
selection limit also leads to miss some fraction of the Blazhko 
variables. 

Concerning the MACHO data, the standard deviations 
of the residuals after the removal of the pulsation and mod- 
ulation components from the light curves are in the 0.03- 
0.20 mag range. This very high noise level makes the de- 
tection of small-amplitude modulations very unlikely. Using 
the light-curve solution of RR Gem and the timings of the 
MACHO observations we generated artificial time series to 
test the detection probability of RR Gem-like modulation 
in the MACHO data. To mimic the noise properties of the 
MACHO data, Gaussian noise was added with 0.03, 0.06, 
and 0.09 mag rms statistics. Fourier analysis of the test data 
(prewhitened for the pulsation) showed that the higher mod- 
ulation peak in the vicinity of the main pulsation frequency 
had an amplitude larger than 2a of the spectrum only if 
the smallest (0.03 mag) rms noise was added. If the rms was 
0.06 mag, there was none among 500 test time series that 
had a modulation peak reaching 2a amplitude. For compar- 
ison, only 5% of the MACHO light curves have rms residual 
smaller than 0.06 mag. 

Our survey is confined only to short-pulsation-period 
variables, which tends to bias the sample in favour of metal- 
rich disk RR Lyraes. Therefore our conclusion that the true 
incidence rate of the modulation is significantly larger than 
it was estimated in the previous surveys is valid only for 
this population of RRab stars. As the MACHO and OGLE 
surveys do not show, however, large differences in the occur- 
rence rates of the modulation at different pulsation periods, 
it is a sound guess that the situation is the same for longer 
period variables, too. To confirm this supposition accurate 
photometric data are needed for longer period RRab stars, 
as well. 



5 THE PHOTOMETRIC METALLICITY OF 
BLAZHKO VARIABLES 

The metallicity of RRab stars can be derived from the pe- 
riod of the pulsation and the tp^i epoch independent phase 
difference of the V light curve with similar accuracy to that 
determined from low di s persion spect roscopic observations 
(IJurcsik fc KovacE|[l996l :lK ovac 3 l2005l ). However, we do not 
know whether this relation also holds for the mean light 
curves of Blazhko variables. If the phase coverage of both the 
pulsation and the modulation cycles is not complete enough 
then the mean light curves of large-modulation-amplitude 
Blazhko variables can be seriously distorted, giving rise to 
spurious value of the Fourier parame t ers. F or instance, us- 
ing the ASAS data of SS For lKovacj (|2005l ) found that the 
photometric metallicity deviates significantly from the spec- 
trosc opic value, however, th e recent, extended photometric 
data (|Kolenberg et al]|2008l ) define a mean light curve with 
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Table 5. Comparison of spectroscopic and photometric [Fe/H] values of Blazhko stars. Variables with small and large modulation 
amplitudes are typeset in italics and boldface, respectively. 



star 


No. of V data 


[Fe/H]sp (weij 


jht)" 


[Fc/H],p " 


[Fc/H]phot 






Suntzeff et al. ( 1994 ~) 


Lavden (1994) 






33 Cnc 


1400 


-0.16 (4.0) 


0.13 (2.5) 


-0.03 


-0.12 


RR Gem 


3000 


-0.14 (4.5) 


-0.13 (2.0) 


-0.14 


-0.16 


DM Cyg 


3000 


-0.16 (1.0) 


0.07 (2.0) 


-0.01 


-0.01 


RY Com 


1000 


-0.64 (0.5) 


-1.38 (1.0) 


-1.13 


-1.16; -1.16; -1.08; -1.23; -l.ie'^ 


AQ Lyr 


1400 


-0.19 (2.0) 






-0.05 


XY And 


2600 




-0.68 (1.0) 




-0.68 


CZ Lac 


8000 


-0.13 (2.0) 






-0.26 


FK Vul'* 


500 




-0.71 (1.0) 




-0.37 


BD Her<* 


400 


-0.13 (2.0) 






-1.05 



" Spectroscopic data transforme d to the metallicity scale of the photometric [Fe/H] and their weights. 
Weighted means as specified in I Jurcsik &: Kovacsl (Il996f) . 

Derived from the mean light curves of the 2007 and 2008, 2009 data before and after JD 2454530 and JD 2454905, respectively. 
Scarced data. 



Fourier parameters well matching the [Fe/H] (P, 1/331) rela- 
tion. 

The extended data we obtained for Blazhko variables 
enable us to check the accuracy of the photometric metallic- 
ities derived from the mean light curves of these variables. 
Among the 14 new Blazhko variables nine have spectroscopic 
abundance determinations. Table [5] lists these variables and 
summarizes their spectroscopic and photometric metallici- 
ties. In this sample there are four and five stars showing 
small- and large-amplitude modulations, respectively. The 
photometric metallicities of the small-modulation-amplitude 
variables calculated from well defined mean fight curves are 
in very good agreement with the spectroscopic values. We 
thus conclude that the photometric metallicity of small- 
modulation- amplitude Blazhko variables is correct, and it 
has the same accuracy as for unmodulated RRab stars if 
the Ught variation is sampled properly. 

The situation for Blazhko variables showing large- 
amplitude modulations is, however, not so obvious. Unfor- 
tunately, the spectroscopic information on these stars is less 
certain than for the small-modulation-amplitude variables. 
For those variables that have good data coverage (AQ Lyr, 
XY And and CZ Lac) the spectroscopic and photometric 
metallicit ies agre e reaso nably well. A similar result was de- 
rived by ISmoled (|2005h too, using published photometric 
and spectroscopic data of five large-modulation-amplitude 
Blazhko stars. For the two large-modulation-amplitude vari- 
ables in our sample, which have only 400-600 measurements 
(FK Vul and BD Her), however, the discrepancy between 
the photometric and spectroscopic metallicities is large. The 
mean Ught curves used to calculate the [Fe/H] phot for FK Vul 
and BD Her are also drawn in the middle panels of Fig. [6] It 
can bee seen that the descending branch of the light curve 
of FK Vul is undersampled, which may result in improper 
shape of the mean light curve. Also, the mean light curve of 
BD Her, especially around minimum light, has an unusual 
shape most probably due to unequal data sampling. It is an 
open question, however, that if the light curves of FK Vul 
and BD Her were better sampled, then their photometric 
[Fe/H] would or would not match their spectroscopic val- 
ues better. The good agreement between the photometric 



and spectroscopic metallicities of those large-modulation- 
amplitude Blazhko variables which have good light-curve 
coverage indicates, however, that if the light variation is well 
sampled then the Fourier parameters of the mean light curve 
can be safely used to calculate the metallicity. Nevertheless, 
to make this statement better grounded, both extended pho- 
tometric data of other large-modulation-amplitude variables 
and spectroscopic metallicities of these variables (e.g. for 
MW Lyr) are needed. 



6 SUMMARY 

• We have shown that with the inclusion of small- 
amplitude modulations the incidence rate of the modula- 
tion increases to about 50% for short-period RRab stars. 
The Konkoly Blazhko Survey showed that, in fact, small- 
amplitude modulations are equally frequent as large ampli- 
tude ones. The distribution of the amplitudes of the modula- 
tion show an increase towards smaller amplitudes according 
to our data. It suggests that the light curves of those RRab 
stars that are now supposed to be stable may also be mod- 
ulated but with very small amplitude, below our present 
detection limit. Moreover, we may have missed detecting 
small-amplitude modulations with long modulation periods 
in our sample. Therefore, it is likely that the incidence rate 
of the modulation we detected is still a lower limit to the 
true fraction of the RRab stars showing the Blazhko effect. 
Modulations with amplitudes of the maximum brightness 
variation in the millimagnitude regime may still have es- 
caped detection. The observations of recent space missions 
(MOST, CoRoT, Kepler) will clear this question soon. 

• There is another embarrassing feature of the Blazhko 
modulation, namely its temporal occurrence as documented 
recentl y in the study of the long term behaviour of RR Gemi- 
norum (|S6dor et al.ll2007bl ') and in the case of RY Com in the 
present paper. The modulation of RR Lyrae i s also know n 
to cease for a short period in every fourth year (|Szeidllll988l ). 
Based on these results there is a good chance that the modu- 
lation is a common feature of the pulsation of fundamental- 
mode RR Lyrae stars, it may be even a universal property 
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of these variables. This possibility warns that understanding 
the Blazhko modulation is crucial in the knowledge of the 
pulsation of RR Lyrae stars. 

• Detecting small-amplitude light-curve modulation of 
RRa b stars i s also relevant from the theoretical point of 
view. iKovacd (|l995l ) used the lack of small-amplitude mod- 
ulations as an argument against the obli que magnetic rota - 
tor model of the Blazhko phenomenon ((Shibahashill200d ). 
He says "Assuming random distribution [of the aspect an- 
gle], we should see some small- and some large-amplitude 
modulations. This is apparently not the case." The discovery 
of small- amplitude modulation contradicts this argument. 
Whatever the cause of the Blazhko modulation, it manifests 
itself in the full possible range of amplitudes. At the same 
time, simultaneous multi-periodic modulation with similar 
modulation amplitudes were also detected in the KBS. The 
significant difference between the modulation periods ob- 
served in these stars is a strong argument against those 
explanations of the Blazhko phenomenon that connect the 
modulation period to the rotation of the star. 

• Based on our photometric data and spectroscopic in- 
formation on Blazhko stars we have found that the photo- 
metric metallicity determined from the mean light curves of 
Blazhko variables are correct if the light curve is properly 
sampled, i.e., the mean light curve is not distorted by uneven 
data sampling of the different Blazhko phases. 

In the Konkoly Blazhko Survey we obtained extended 
data on 30 RRab stars with fundamental periods shorter 
than 0.5 days, 14 of them showing the Blazhko effect. This is 
a significantly larger incidence rate of the modulation than 
found previously. We cannot exclude the possibility, how- 
ever, that the large percentage of Blazhko stars in the sample 
of the KBS is related to the short pulsation period of these 
stars. If this were the case, the incidence rate of the modu- 
lation would also depend on the metallicity of the stars as 
short period RRab stars are, in general, less metal poor than 
longer period variables. In order to decide whether there is a 
period and/or metallicity dependence of the occurrence rate 
of the modulation we initiated a similar survey focusing on 
RRab stars of longer periods in 2009. 
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